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Abstract 

The  modem  world  has  an  inereasing  demand  for  wireless  multiple  access 
communications;  direct-sequence  spread-spectmm  multiple  access  (DS/SSMA)  systems 
comprise  many  of  these  communication  systems.  A  better  understanding  of  multiple 
access  interference  (MAI)  effects  on  DS/SSMA  system  performance,  specifically  their 
impact  on  overall  system  bit  error  rate  (BER),  enables  system  designers  to  minimize  MAI 
degradation  and  produce  greater  DS/SSMA  system  capacity. 

This  research  characterizes  MAI  effects  on  DS/SSMA  system  performance 
through  simulation  in  Matlab®,  and  explores  the  impact  of  multiple  access  code  selection, 
chip  waveform  shaping,  and  multiple  access  code  length  on  BER  for  both  synchronous 
and  asynchronous  multiple  access  networks.  In  addition,  the  simulated  DS/SSMA  model 
permits  rapid  research  into  the  effects  of  additional  factors  on  BER. 

Prior  to  experimental  testing,  model  validation  is  conducted  through  single  user 
trials  and  by  comparison  with  existing  research  for  similar  system  designs.  Eor 
synchronous  and  asynchronous  networks.  Gold  coding  improves  BER  by  7.5  and  4.0  dB, 
respectively,  relative  to  aperiodic  random  spreading  codes.  Synchronous  network  results 
show  that  chip  waveform  shaping  provides  no  significant  BER  improvement  for  the 
Blackman  or  Eanczos  shapes.  However,  asynchronous  network  results  show  a  potential 
BER  improvement  for  Blackman  and  Eanczos  shapes.  Increasing  code  length  from  31  to 
511  resulted  in  a  7.5  dB  BER  improvement.  Collectively,  these  results  directly  relate 
changes  in  BER  to  waveform  cross-correlation  statistics. 


XI 


MULTIPLE  ACCESS  INTEREERENCE  CHARACTERIZATION  EOR 


DIRECT-SEQUENCE  SPREAD-SPECTRUM  COMMUNICATIONS 
USING  CHIP  WAVEEORM  SHAPING 

I.  Introduction 

1.1  Motivation  for  Multiple  Access  Interference  Characterization 

The  ever  increasing  demand  for  world-wide  multiple  access  wireless 
communications  drives  the  need  to  maximize  the  current  system  capability  and 
transmission  capacity.  Many  current  systems  employ  direct-sequence  spread-spectrum 
(DSSS)  techniques  to  enable  multiple  access  capability.  A  clearer  understanding  of  how 
multiple  users  effect  overall  system  performance,  through  characterization  of  the  multiple 
access  interference  (MAI),  enables  more  efficient  use  of  current  systems  and  better 
designs  for  future  systems. 

Although  a  large  body  of  research  exists  on  how  MAI  impacts  direct-sequence 
spread-spectrum  multiple  access  (DS/SSMA)  performance,  most  of  this  work  relies  on 
analytical  approximations.  Additionally,  there  are  numerous  factors  that  can  affect  MAI 
contributions  in  a  system  between  workload  and  environment  (e.g.,  number  of 
simultaneous  transmitters,  type  of  multiple  access  coding,  and  code  length).  Current 
approximations  only  account  for  a  limited  number  of  these  factors  and  can  require 
extensive  recalculation  when  factors  are  changed  or  added.  The  state  of  existing  research 
into  MAI  in  DS/SSMA  systems  leaves  open  a  need  for  a  representative  system  model  that 
is  easily  modified  to  account  for  different  factors. 
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This  need  lends  itself  to  eharaeterization  of  MAI  in  DS/SSMA  systems  through 
simulation.  Simulation  of  DS/SSMA  system  performanee  provides  two  main  benefits. 
First,  simulation  allows  for  rapid  testing  of  the  effeets  that  numerous  faetors  ean  have  on 
system  performanee  and  sueh  faetors  ean  be  ehanged  or  others  added  relatively  easily. 
Seeond,  simulation  enables  verifieation  of  future  approximations  through  a  vehiele  that 
more  elosely  represents  the  aetual,  physieal  eommunieation  system. 

1.2  Problem  Statement  and  Scope 

To  inerease  the  eapaeity  and  eapability  of  DS/SSMA  systems,  a  greater 
understanding  of  MAI  effeets  (as  a  funetion  of  multiple  faetors)  on  system  performanee 
is  required.  Simulation  of  a  DS/SSMA  system  model  provides  the  ability  to  eharaeterize 
the  effeets  of  various  design  faetors  and  permits  verifieation  of  existing  and  future 
analytieal  approximations.  This  work  provides  modeling  and  simulation  results  for  a 
representative  DS/SSMA  system  and  inereases  the  understanding  of  partieular  faetors' 
impaet  on  MAI.  Speoifieally,  this  work  simulates  the  effeet  of  multiple  aeeess  eode 
seleetion,  ehip  waveform  shaping,  and  multiple  aeeess  eode  length  on  bit  error  rate 
(BER)  in  a  DS/SSMA  system.  Furthermore,  this  researeh  provides  a  simulation  base 
eapable  of  relatively  easy  modifieation  to  evaluate  the  impaet  of  additional  faetors  on 
BER. 

1 . 3  Methodology 

This  researeh  simulates  the  transmitter,  ehannel,  and  reeeiver  of  a  DS/SSMA 
system  using  Matlab®.  The  simulation  eontains  all  elements  of  random  binary  data 
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generation,  data  modulation,  multiple  aceess  coding,  transmission,  reception, 
despreading,  and  communication  symbol  detection  and  estimation.  All  environmental 
and  system  level  factors  that  potentially  impact  BER  are  contained  within  the  simulation 
and  easily  modified  to  support  additional  research.  Theoretical  performance  models  and 
previous  research  provide  the  basis  for  system  model  verification.  Simulation  is  used  to 
estimate  the  effect  of  multiple  access  code  selection,  chip  waveform  shaping,  and 
multiple  access  code  length  on  BER. 

1.4  Thesis  Organization 

This  thesis  contains  five  chapters.  Chapter  1  provides  an  introduction  on  the  need 
for  MAI  characterization  for  DS/SSMA  systems.  Chapter  2  outlines  relevant  background 
information  on  existing  BER  approximations  and  the  factors  to  be  explored  in  research. 
Research  methodology  is  outlined  in  Chapter  3.  Verification  and  experimental  testing 
results  are  provided  in  Chapter  4.  Chapter  5  contains  a  summary  of  contributions  and 
findings  and  outlines  possible  future  research.  One  appendix  is  provided  containing  the 
Matlab®  code  associated  with  system  simulation. 


1-3 


II.  Background 


2. 1  Direct  Sequence  Spread  Spectrum  Communications 

Modern  digital  communication  techniques  have  kindled  widespread  appeal  and 
demand  for  improving  existing  and  future  communications  systems.  This  increased 
demand  has  forced  regulating  organizations  to  establish  limits  on  what  portions  of  the 
spectrum  specific  applications  can  use.  The  result  is  limited  bandwidth  for  all 
communications  applications.  The  desire  to  get  more  use  out  of  the  available  bandwidth 
has  driven  various  techniques  to  increase  the  number  of  users  or  applications 
simultaneously  occupying  a  specific  spectral  region.  One  such  technique  is  direct- 
sequence,  spread-spectrum  multiple  access  (DS/SSMA). 

2. 1. 1  Multiple  Access  Interference. 

Direct-sequence  spread-spectrum  (DSSS)  techniques  are  commonly  used  to 
implement  multiple  access  communications.  Existing  systems  using  DSSS  techniques 
include  the  Global  Positioning  System  (GPS)  and  IS-95  digital  cellular  phone  system 
[Peterson,  Ziemer,  and  Borth,  1995].  DSSS  coding  modulates  each  communication 
symbol  with  a  waveform  consisting  of  multiple  chip  intervals.  The  chip  modulating 
waveform  has  a  particular  pulse  shape  with  either  a  positive  or  negative  orientation  as 
determined  by  the  DSSS  code.  The  undesired  signals  make  up  the  multiple  access 
interference  (MAI)  that  increases  with  an  increasing  number  of  users.  For  perfectly 
orthogonal  DSSS  coding  the  MAI  term  is  identically  zero  and  the  system  functions  as  if 
there  were  only  a  single  user  present,  i.e.,  the  desired  user.  Although  the  orthogonal 
coding  is  optimal,  the  number  of  required  users,  mathematical  code  limitations,  the 
computational  complexity  required  to  generate  large  numbers  of  orthogonal  codes,  and 
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the  desire  to  analytieally  model  systems  typieally  neeessitates  the  selection  of  non- 
orthogonal  codes  in  either  the  design  or  analysis  of  DS/SSMA  systems  [Gerantiotis  and 
Ghaffari,  1991].  DS/SSMA  systems  rely  on  code  cross-correlation  characteristics  to 
reject  undesired  signals  and  minimize  MAI.  The  ability  to  reject  undesired  signals  is 
known  as  multiple  access  interference  suppression.  By  reducing  MAI  in  a  given  system, 
more  users  can  simultaneously  access  the  same  link. 

2.1.2  Signal  and  System  Model. 

The  works  of  Yao  and  Pursley  present  a  standard  DS/SSMA  system  design.  The 
system  consists  of  K  users  with  each  using  BPSK  modulation  for  their  respective  data. 
The  data  modulated  waveform  bi(t)  for  each  user  can  be  represented  by  [Yao,  1977] 

CO 

biid)=  -  nTJ  (2.1) 

n=—cc 

where  bi„  is  a  sequence  of  elements,  b^n  s  [-1,1],  and  Pt(*)  is  a  unit  height  rectangular 
pulse  of  duration  Ts,  the  symbol  duration.  Spreading  code  modulation  ak(t)  is  defined  as 
[Pursley,  1977] 

00 

«yt(0=  Pyit  -  jT^)  (2.2) 

j=-«> 

where  is  a  sequence  of  elements  for  the  user,  a,-  8  [-1,1],  each  chip  interval  has 
duration  J/,  and  the  spreading  code  period  is  N»Tc.  The  most  common  pulse  shape  used 
over  each  chip  interval  is  rectangular,  as  represented  by  (t-  jT^  )  in  (2.2).  One  of 

the  most  common  spreading  techniques  involves  multiplying  each  communication 
symbol  by  N  coded  replications  of  the  chip  waveform  over  one  full  code  period  per 
symbol  interval.  This  particular  technique  is  the  basis  for  DSSS  systems  considered  in 
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this  research.  Using  the  data  and  spreading  code  modulation  of  (2.1)  and  (2.2), 
respectively,  the  transmitted  signal  for  the  user  is  defined  as 

H  (0  =  ViP  (O^i  (0  cos(®,^  +  P J  (2.3) 

For  the  equal  power  cases  considered  here,  represents  the  transmitted  power  of  all 
signals,  cos((»^^  +  is  the  phase  modulated  carrier  of  the  user,  cOc  is  the  modulation 

frequency,  and  6k  is  the  phase  delay  of  the  user  [Pursley,  1977]. 

The  total  received  multiple  access  signal  is  the  sum  of  K  transmitted  signals  and 
the  channel  noise  and  may  be  expressed  as 

K 

r{t)  =  ^  VZP  {t  -  )/)^ (^  -  U ) cos{(o^t  +  ^k)  +  «(0  (2.4) 

k=\ 

where  Tk  and  are  the  time  delay  and  received  phase,  respectively,  of  the  user.  The 

function  n(t)  represents  thermal  channel  noise  and  is  assumed  to  be  additive  white 
Gaussian  noise  (AWGN)  having  two-sided  spectral  density  of  No/2.  All  subsequent 
derivations  presented  in  this  work  assume  the  desired  user's  signal  has  zero  time  delay 
and  zero  received  phase.  All  interfering  users  have  time  delays,  relative  to  the  desired 
user,  in  the  interval  [0,  Ty]  and  relative  phases  in  [0,  2n].  In  the  special  case  of  a 
synchronous  network  all  interfering  users  have  zero  time  delay. 

2.2  Analytical  Bit  Error  Rate  (BER)  Approximation 

The  work  of  [Pursley,  1977]  provides  an  expression  for  bit  error  rate  (BER)  in 
terms  of  the  number  of  network  users  {K)  and  code  length  {N).  Using  a  conventional 
matched  filter  design,  the  received  signal  is  despread  and  demodulated  via  correlation 
which  generates  the  test  statistic  Z,  given  as  [Pursley,  1977] 
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(2.5) 


T 

Z;  =  J  r{t)a^  {t)  cos(cL)J)dt . 

0 

Assuming  (»c  >>  Ti^,  the  double  frequency  component  can  be  ignored  and  Z  can  be 
rewritten  as 


k=i 

k*i 


T 

+  J  n{t)a^  {t)  cos((i)j)dt 
0 


(2.6) 


where  Rj^j  (t)  and  Rj^  i  (x)  are  defined  as  continuous-time,  partial  cross-correlation 
functions  and  expressed  as  [Pursley,  1977] 


RkM)  =  \ak{t-'!^)ai{t)dt 
0 
T 

h,ii'^)  =  \ak{t-r)a^{t)dt 


(2.7) 


The  partial  cross-correlations  R^.  ■  (r)  and  7?^^  ■  (r)  of  (7)  account  for  mismatch  between 

symbol  transition  boundaries  of  the  desired  user  and  the  k‘^  user  resulting  from  the 
user’s  time  delay,  and  the  user  is  an  interferer.  R^. .  (r)  correlates  the  end  of  the  first 

symbol  of  the  A:'*  user  that  falls  within  the  desired  user’s  symbol  interval,  and  R,^  ,  (r) 

correlates  the  beginning  of  the  next  subsequent  symbol  of  the  k**'  user  which  completes 
the  rest  of  the  desired  user’s  symbol  interval. 

Rj^l(  X )  and  Rj^i(  x )  can  be  written  in  terms  of  the  discrete  aperiodic  cross¬ 
correlation  function  Ck,i(l)  defined  as  [Pursley,  1977] 
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(2.8) 


Q.,(/)  = 


N-l-l 


0</<  A-1 

J=0 

N-l+l 

(k)  i 
l^aYiUj 

l-A</<0 

j=0 

0, 

|/|>A. 

Substituting  Ckj(l)  of  (2.8)  into  (2.7)  results  in  the  following  alternate  expressions  for 

Rkj(T): 


Rk,A^)  =  -  NK  +  + 1  -  A^)  -  Q,,(/  -  N)]  •  (r  -  ITJ 

hM  =  C,Al)T^  +  [C, ,(/  + 1)  -  Q, ,(/)]  •  (r  -  /T;) 


(2.9) 


To  establish  the  reeeived  signal-to-interferenee-plus-noise  ratio  (SINR),  the 
power  of  the  noise  plus  interfering  signals  is  ealeulated  as  the  varianee  of  Z,  defined  in 
(2.6)  over  one  symbol  interval  and  is  given  by  [Pursley,  1977] 

D  A  K  T 


Var{ZA  = 


4T. 


4T . 


Sjh  k/ (^)  +  {r))dr  +  NJ  /  4 

k  =  \o 
kiti 

K  v-i(^+9T  N 

2  2  I  +  4,,7r))rfr  +  f^J/4. 


(2.10) 


k=l  1=0  IT 
k*i 


In  this  ease,  the  statistieal  expeetation  is  taken  with  respeet  to  variables  (pk,  T^k,  and  bt- 
Given  the  desired  user’s  data  remains  eonstant  over  the  symbol  interval  under 
eonsideration,  (2.10)  only  eontains  eontributions  from  ehannel  noise  and  MAI. 

Substituting  the  7?^^  ,Yr)and  Ri^  i(t)  expressions  from  (2.9)  into  (2.10)  results  in  a 
simplified  varianee  given  by 


Var{ZA  = 


(  pj'^  ^ 


K 

2‘ 

k=\ 

\kH 


kj 


+  NJ/4 


(2.11) 
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where  [Pursley,  1977] 


N-l 


'kj 


=  ^  (  C,/ (I -N)  +  Q,,  (/  -  N)C,,  (I -N  +  \)  +  C,/  (I -N  +  \)  (2.12) 


+  Q/  (/)  +  C, ,  (/)C, ,  (/  +  !)  +  C,/  (/  + 1) 


The  SINR  is  then  determined  by  dividing  ■yJP/2T  by  the  rms  noise  of  ^Var{ 


[Pursley,  1977],  yielding; 


SINR  =  ] 


(6N^r''Zn,+ 


k=l 

k^i 


-1/2 


(2.13) 


To  faeilitate  preliminary  system  design,  it  is  shown  in  [Pursley,  1977]  that  for  random 
spreading  sequenees  (2.13)  simplifies  to 


K-1 

SINR  =  < - + 


-1/2 


3N  2E^ 


(2.14) 


where  the  (K-I)/(3N)  term  results  from  taking  the  expeeted  value  of  the  summation  in 
(2.13)  given  random  sequenees  are  employed.  The  equation  for  probability  of  bit  error 
(Pb)  for  a  BPSK  system  without  interferers  is  [Peterson,  Ziemer,  and  Borth,  1995] 


Pb=Q 


2E, 


N, 


(2.15) 


o  J 


Using  (2.14)  an  approximate  Pb  value  ean  be  ealeulated  for  a  DS/SSMA  system  using 
(2.15)  as 

Ps=Q{SINR).  (2.16) 


It  is  useful  to  note  that  if  only  one  user  is  present  {K  =  1),  (2.14)  identieally  simplifies  to 
the  g-funetion  argument  of  (2.15),  the  eorreet  result  for  the  single  user  ease. 
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2.2. 1  BER  for  Asynchronous  and  Synchronous  Systems. 

Using  characteristic  functions  for  DSSS  signals,  Gerantiotis  and  Ghaffari  have 
verified  the  BER  approximation  results  of  Pursley  [Gerantiotis  and  Ghaffari,  1991; 
Pursley,  1977]  and  redefined  (2.15)  via  derivation  as 

SINR  =  \^^m^+^\  (2.17) 

1  N  ^  2E,\ 

where  =  1/3  for  reetangular  pulses  and  is  defined  as 


(2.18) 


0  0 

In  (2.18)  R^  and  R^  are  partial  autoeorrelation  funetions  of  the  pulse  shape  defined  as 
[Gerantiotis  and  Ghaffari,  1991] 


T 


(2.19) 


This  derivation  elarifies  that  (2.14)  is  specifieally  applieable  to  the  SINR  approximation 
for  asynchronous  signals  using  reetangular  shaped  ehip  waveforms  [Gerantiotis  and 
Ghaffari,  1991].  Continuing  with  this  development  a  BER  approximation  for 
synchronous  DS/SSMA  systems  using  reetangular  shaped  ehip  waveforms  is  found  to  be 
[Gerantiotis  and  Ghaffari,  1991] 

=  +  ^  .  (2.20) 
1  2/V  2E,  I 

It  has  been  shown  that  BER  ean  be  approximated  for  asynchronous  and  synchronous 
systems  by  substituting  (2.14)  and  (2.20),  respeetively,  into  the  g-funetion  argument  of 
(2.16)  [Gerantiotis  and  Ghaffari,  1991]. 
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2.2.2  Accuracy  of  BER  Approximations. 

Using  the  SINR  approximations  of  (2.14)  and  (2.20),  and  their  own  eharaeteristie 
funetion  approximations  for  SINR,  Gerantiotis  and  Ghaffari  ealeulated  the  BER  for  both 
the  synehronous  and  asynehronous  eases  using  both  random  aperiodie  spreading  eodes 
and  deterministie  m-sequenee  codes  [Gerantiotis  and  Ghaffari,  1991].  For  all  cases 
considered,  parameter  values  of  if  =  3  and  A^=  3 1  were  used  with  rectangular-shaped 
chip  waveforms.  Their  results  for  PSK  modulation  are  shown  in  Table  2. 1 .  All  values 
with  superscript  G  are  approximations  resulting  from  (2.14)  and  (2.20),  and  those  without 
a  superscript  result  from  Gerantiotis  and  Ghafferi's  characteristic  function  approximation. 


Table  2.1.  Error  Probabilities  for  DS/SSMA  Systems  (if  =  3,  A  =  31) 
[Gerantiotis  and  Ghaffari,  1991] 


Random  Sequences 

m-Sequences 

Synchronous 

Asynchronous 

Synchronous 

Asynchronous 

N.. 

Pf 

P 

P 

K 

P. 

p: 

P. 

8 

1.37x10'^ 

1.57x10'^ 

8.15x10'“ 

9.13x10'“ 

1.03x10'® 

8.04x10'“ 

8.13x10'“ 

8.49x10'“ 

10 

2.45x1 0'"' 

3.71x10''* 

9.21x10'® 

1.44x10' 

1.43x10'“ 

5.99x10® 

9.17x10'® 

1.08x10'“ 

12 

3.77x10'^ 

9.73x10'® 

7.07x10'® 

2.55x10'® 

1.54x10'® 

1.28x10'® 

7.02x10'® 

1.08x10'® 

14 

5.98x10'® 

3.19x10'® 

4.46x10'’ 

7.69x10'® 

153x10'** 

3.96x10'" 

4.40x10'’ 

1.05x10'® 

16 

1.16x10'® 

1.35x10'® 

3.01x10'** 

44.1x10'’ 

1.78x10'’ 

5.98x10'*’ 

2.95x10'** 

1.06x10'’ 

Four  observations  from  Table  2.E  are  important  for  this  research.  First,  the 
analytical  BER  approximations  based  on  (2.14)  and  (2.20)  were  found  to  be  optimistic 
for  random  codes  in  the  synchronous  and  asynchronous  cases,  and  for  deterministic  codes 
in  the  asynchronous  case.  However,  the  BER  approximations  based  on  (2.14)  and  (2.20) 
are  very  conservative  for  synchronous  systems  with  deterministic  spreading  codes. 
Second,  for  random  sequences  the  asynchronous  case  outperforms  the  synchronous  case 
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due  to  the  inereased  randomization  and  averaging  provided  by  the  asynehronous  ease 
[Gerantiotis  and  Ghaffari,  1991].  Third,  for  deterministic  sequences  the  synchronous 
case  outperforms  the  asynchronous  case,  reflecting  the  fact  that  deterministic  codes  are 
typically  designed  to  maximize  the  ratio  of  autocorrelation  to  cross-correlation  for  the 
synchronous  case  [Peterson,  Ziemer,  and  Borth,  1995].  Finally,  in  both  synchronous  and 
asynchronous  cases  the  deterministic  codes  outperformed  the  random  codes.  This  result  is 
more  pronounced  in  the  synchronous  case  [Gerantiotis  and  Ghaffari,  1991]. 


2.3  Chip  Waveform  Shaping 

To  reduce  MAI  effects  in  DS/SSMA  systems,  the  impact  of  using  chip  waveform 
shapes  other  than  rectangular  is  explored.  This  approach  to  reducing  MAI  is  based  on 
decreasing  the  cross-correlation  statistics  between  spreading  codes  as  a  result  of  variation 
in  the  chip  waveform  shape.  The  pulse  shapes  considered  in  this  research  include 
rectangular,  half-sine,  raised-cosine,  Blackman,  Kaiser,  and  Lanczos.  As  defined  in  [Kok 
and  Do,  1997]  each  of  these  pulse  shapes  can  be  analytically  represented  as  follows: 
a)  Blackman 


y/{t)  =  k 


f 

0.42 -0.5  cos 

u(t) 


(2.21) 


b)  Half-sine 


(//(t)  =  V2  sin 


M 


yT.j 


u{t) 


(2.22) 
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c)  Kaiser 


d)  Lanczos 


e)  Raised-cosine 


=  k 


t-TJl 


u{t) 


\j/{t)  -  k 


smilTU) 

Ijzt 


u(t) 


( 

1-cos 

V 


ItU' 


u{t) 


(2.23) 


(2.24) 


(2.25) 


where  Tc  is  the  chip  duration,  u(t)  is  the  unit  step  function,  Iq  is  a  zero-order  modified 
Bessel  function,  and  A:  is  a  real-valued  scaling  constant. 

2.3.1  Random  Coding  with  Chip  Waveform  Shaping 

Using  the  system  model  from  (2. 1-2.4)  and  chip  waveform  definitions  of  (2.21), 
(2.22),  and  (2.25),  Lehnert  and  Cho  investigated  the  impact  of  waveform  shaping  on 
DS/SSMA  system  performance  using  aperiodic  random  spreading  sequences  [Lehnert, 
2002;  Cho  and  Lehnert,  1999].  All  results  are  for  systems  with K=9  total  users  and  A  = 
3 1  for  the  processing  gain.  Each  user  has  equal  energy,  and  the  time  delays  of  the  8 
interfering  users  are  normalized  to  one  symbol  interval  relative  to  the  desired  user  and 
modeled  as  uniformly  distributed  random  variables. 

The  average  BERs  are  computed  using  (2.14),  where  a  conditional  Gaussian 
approximation  (CGA)  is  used  to  estimate  the  interference  term  when  estimating  SINR, 
which  is  the  argument  of  (2.16).  The  CGA  is  used  here  in  place  of  standard  Gaussian 
approximation  instrumental  in  the  derivation  of  (2.14)  and  (2.20)  because  the  CGA  more 
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accurately  approximates  the  MAI  effeet  on  eommunieation  performanee  in  a  DS/SSMA 
system  [Cho  and  Lehnert,  1999].  Based  on  the  CGA  assumption,  the  Blaekman 
waveform  shape  was  predieted  to  provide  the  best  performanee  relative  to  the  half-sine, 
raised-eosine,  and  reetangular  shapes.  Results  of  this  estimation  proeess  are  reprodueed 
in  Fig.  2. 1 .  for  half-sine,  raised-eosine,  and  Blaekman  pulse  shapes  at  EbINo  values 
ranging  from  0  to  30  dB.  The  data  elearly  illustrates  the  Blaekman  ehip  waveform  shape 
outperforms  all  other  waveform  shapes  eonsidered. 


Iff' 


10'^ 

Average 
BER 

Iff^ 


Iff"' 


0  5  10  15  20  25  30 

Eb/No  (dB) 

Figure  2.1.  Random  Coding  with  Chip  Waveform  Shaping;  Average  BER  for  half- 
sine,  raised-eosine  and  Blaekman  shapes  [Cho  and  Lehnert,  1999]. 

2.3.2  Gold-Coding  with  Chip  Waveform  Shaping. 

The  work  of  Kok  and  Do  analytically  explores  the  impact  of  using  chip  waveform 
shapes  on  DS/SSMA  BER  performanee.  All  shapes  evaluated  are  time-limited  and 
normalized  for  the  analysis.  As  in  the  random  eoding  ease,  the  ehip  waveform  shapes 
evaluated  were  reetangular,  Blaekman,  half-sine,  Kaiser,  Lanezos,  and  raised-eosine. 
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Using  Gold  codes  of  length  127  and  51 1,  while  inereasing  number  of  users  from  5  to  70, 
the  SINR  is  ealeulated  for  eaeh  pulse  shape  case  based  on  aperiodie  autoeorrelation 
funetions  of  (2.8)  [Kok  and  Do,  1997].  This  SINR  value  is  then  used  as  the  argument  of 
(2.16)  to  estimate  BER  for  each  case,  where  the  g-funetion  is  evaluated  with  the 
algorithm  deseribed  in  [Pari,  1980].  Results  show  that  all  non-reetangular  ehip  waveform 
shapes  outperform  the  reetangular  ease.  Speeifieally,  the  system  employing  the  Lanezos 
shape  aehieves  the  best  BER,  followed  in  order  by  Kaiser,  Blaekman,  raised-eosine,  half¬ 
sine,  and  reetangular  [Kok  and  Do,  1997]. 

2.4  Summary 

The  goal  of  this  researeh  is  to  eharaeterize  the  effeets  of  pulse  shaping  on 
DS/SSMA  through  simulation.  This  ehapter  provides  a  deseription  of  the  DSSS  system 
design  and  the  MAI  term  eaused  by  multiple  simultaneous  transmitters.  An 
approximation  for  expeeted  BER  for  both  synehronous  and  asynehronous  systems  is 
derived.  Diseussion  of  the  aeeuraey  of  this  BER  approximation  enables  a  more  informed 
eomparison  to  the  results  of  this  researeh.  Einally,  the  basies  of  pulse  shaping  as  a  means 
of  BER  improvement  is  diseussed  for  random  and  Gold  eodes,  and  examples  of 
approximated  performanee  improvement  are  provided  as  a  baseline  for  eomparison  to 
this  researeh. 
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III.  Methodology 


3.1  Problem  Definition 

This  research  characterizes  the  multiple  access  interference  (MAI)  for  direct- 
sequence,  spread-spectrum  multiple  access  (DS/SSMA)  systems  through  simulation,  and 
investigates  the  impact  of  code  selection,  pulse  shaping,  and  code  length  on  bit  error  rate 
(BER).  Multiple  access  capabilities  are  measured  in  terms  of  BER  for  given  signal-to- 
noise  ratios  as  related  to  the  average  energy  per  bit  of  the  desired  user  {Eb)  divided  by  the 
background  noise  power  spectral  density  {No).  The  effects  that  parametric  changes  have 
on  multiple  access  BER  are  also  considered.  Previous  analytical  work  characterizes  BER 
in  terms  of  Et/No,  the  number  of  users,  and  the  ability  of  pulse  shaping  over  a  chip 
interval  to  reduce  MAI  levels.  This  research  develops  and  uses  a  Matlab®  simulation  to 
model  DS/SSMA  system  performance  and  to  investigate  the  impact  code  selection  and 
pulse  shaping  has  on  BER.  Additionally,  the  simulation  provides  a  means  for  predicting 
the  impact  of  increasing  code  length  in  systems  employing  random  and  Gold  codes  in 
conjunction  with  the  pulse  shaping.  Simulation  results  are  verified  against  the  single  user 
communications  baseline  and  research  results  outlined  in  Chapter  2.  Once  verified,  the 
code  is  used  to  characterize  the  impact  of  synchronization,  code  selection,  code  length, 
and  pulse  shaping  on  BER  for  a  DS/SSMA  system. 

3.2  Approach 

Simulation  results  are  reported  as  BER  vs.  Eb/No  curves  for  both  synchronous  and 
asynchronous  systems  with  multiple  users,  using  random  binary  and  Gold  coding 
techniques,  and  implementing  chip  waveform  shaping.  Each  system  user  transmits 
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communication  waveforms  using  binary  phase  shift  keyed  (BPSK)  baseband  modulation 
and  coded  spreading  waveforms  consisting  of  the  desired  pulse  shape  over  the  chip 
intervals.  All  signals  are  assumed  to  be  reeeived  with  equal  power  and  additive  white 
Gaussian  noise  (AWGN)  power  levels  are  adjusted  to  achieve  the  desired  Eh/No  value. 
Each  chip  interval  has  duration  Tc  and  spreading  eodes  are  periodie  with  length  N.  Thus, 
a  colleetion  of  coded  ehip  waveforms  eomprises  one  symbol  duration  Ts  and  equals  NTc. 
Simulation  validation  consists  of  two  stages,  ineluding  1)  single  user  results  are  generated 
and  compared  with  (2-15)  and  2)  results  are  compared  to  the  analytical  findings  of 
Gerantiotis  and  Ghaffari  [Gerantiotis  and  Ghaffari,  1991].  Following  code  validation, 
simulation  results  are  used  to  explore  the  impaet  of  pulse  shaping,  code  seleetion,  and 
code  length  on  BER  for  DS/SSMA  systems. 

3.3  System  Boundaries 

The  system  of  interest  is  a  direet-sequenee  spread-spectrum  (DSSS)  receiver 
using  BPSK  modulation  (data  and  spreading)  for  multiple  access  communications. 

Figure  3.1  shows  a  diagram  of  the  physical  system  whieh  eonsists  of  a  reeeive  antenna,  a 
radio  frequeney  (RE)  bandpass  filter,  a  despreading  mixer,  an  intermediate  frequeney  (IF) 
bandpass  filter,  and  the  data  phase  demodulator.  This  research  simulates  DSSS  receiver 
performance  in  a  multiple  access  environment  using  Matlab®.  Each  user  transmits 
eommunieation  waveforms  consisting  of  1)  BPSK  data  modulation  and  2)  BPSK 
spreading  modulation  having  the  desired  pulse  shape  over  each  chip  interval.  All  signals 
are  reeeived  with  equal  power  and  AWGN  power  levels  are  adjusted  to  aehieve  the 
desired  Eh/No  value.  Each  chip  interval  has  duration  Tc  and  spreading  codes  are  periodic 


3-2 


with  length  N  such  that  the  symbol  duration  equals  NTc.  The  ability  of  the 
despreading  mixer,  IF  filter  and  demodulator  to  correctly  estimate  communication 
symbols  (which  are  subsequently  mapped  to  bits)  from  the  desired  user,  in  the  presence 
of  other  system  users,  is  the  focus  of  this  research. 

V 


Despreading 

Mixer 


RF 

IF 

Data  Phase 

Filter 

Filter 

Demodulator 

Figure  3.1.  BPSK  DSSS  Receiver  [Peterson,  Ziemer,  and  Borth,  1995]. 


3.4  System  Services 

DSSS  receivers  are  designed  to  receive  signals  and  estimate  bits  sent  from  the 
desired  user  while  rejecting  interfering  signals  of  all  undesired  users.  By  doing  so,  DSSS 
receivers  enable  multiple  users  to  simultaneously  communicate  within  the  same  spectral 
region.  Over  each  received  symbol  interval,  the  receiver  estimate  results  in  one  of  two 
possible  outcomes,  either  1)  a  symbol  (bit)  is  received  and  estimated  correctly  (correct 
condition),  or  2)  a  symbol  (bit)  is  received  and  estimated  incorrectly  (error  condition). 

The  source  of  making  errors  per  outcome  2  can  be  either  background  noise  or 
interference  from  multiple  users  (MAI).  Although  AWGN  background  noise  is  present 
for  all  simulations  and  analyses  conducted  in  this  research,  the  primary  focus  of  this  work 
is  on  errors  due  to  MAI. 
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5.5  Performance  Metrics 

The  introduction  of  errors  in  signal  detection  and  estimation  for  signals  received 
through  a  physical  medium  is  inevitable  and  these  errors  are  even  more  likely  to  occur 
with  MAI  present.  The  percentage  of  errors  occurring  in  DSSS  receiver  processing  is 
quantified  using  bit  error  rate  (BER  or  Pb),  defined  as  the  number  of  demodulated  bits 
that  are  incorrectly  estimated  at  the  receiver  divided  by  the  total  number  of  bits 
transmitted.  The  required  system  BER  establishes  a  lower  bound  on  the  percentage  of 
errors  a  receiver  can  make  and  still  correctly  function.  This  research  uses  BER  as  the 
primary  metric  to  measure  system  performance.  Eundamental  to  BER  determination  is 
the  ratio  of  available  signal  energy  Eh  to  total  interfering  energy.  In  this  case,  the  total 
interfering  energy  consists  of  channel  noise  energy,  as  established  by  No,  and  interference 
due  to  multiple  users  (A/).  The  resultant  energy-to-noise  ratio  can  be  expressed  as 


Eh 

No+N/ 


(3.1) 


Eor  a  system  using  BPSK  data  modulation,  BER  is  directly  related  to  this  ratio  and  is 


where  Q  is  the  complimentary  error  function  and  Nt  =  No  +  Nj  [Eehnert,  2002].  The 
existing  analytical  approximations  estimate  Et/Nr  values  which  are  converted  to  BER  for 
comparison  with  the  simulated  results  using  (3.2). 
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3.6  Parameters 


3.6.1  System 

System  parameters  whieh  affeet  the  demodulator’s  ability  to  eorreetly  estimate 
symbols  from  the  desired  user  inelude: 

•  Channel  propagation  oharaeteristies 

•  Reeeiver  antenna  gain  (or  loss) 

•  Speetral  width  and  eenter  frequencies  of  the  RF  and  IF  filters 

•  User  of  interest  (one  of  total  system  users) 

•  Type  of  spreading  code  used 

These  system  parameters  reflect  the  physical  receiver  as  simulated.  As  indicated  earlier, 
the  channel  is  assumed  to  be  an  ideal  AWGN  channel  with  a  two-sided  constant  power 
spectral  density  of  No/2.  Receive  antenna  gain  affects  the  total  amount  of  received  power 
(signal,  noise,  and  interference).  For  this  work,  it  is  assumed  that  received  antenna  gain 
affects  all  signal,  noise,  and  interference  terms  equally  and  is  therefore  ignored.  The 
spectral  width  and  center  frequencies  of  the  RF  and  IF  bandpass  filters  impact  the  amount 
of  signal  power  processed  by  the  system.  Filter  center  frequency  is  determined  by  which 
user  is  designated  as  the  desired  user,  and  the  fdter’s  spectral  width  and  shape  are 
determined  by  the  spreading  code  employed.  For  all  simulations,  it  is  assumed  that  the 
center  frequency  and  spreading  code  are  perfectly  selected  for  the  desired  user.  The  type 
of  spreading  codes  employed  impact  the  level  of  MAI  due  to  varying  autocorrelation  and 
cross-correlation  characteristics  of  different  codes.  Random  codes  are  employed  for  both 
single  user  and  multiple  access  validation  tests.  This  is  appropriate  for  the  single  user 
validation  because  single  user  performance  is  independent  of  code  selection.  Random 
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codes  are  employed  for  multiple  aeeess  validation  to  permit  direet  eomparison  with  the 
system  setup  used  by  Gerantiotis  and  Ghaffari  [Gerantiotis  and  Ghaffari,  1991].  Both 
random  and  Gold  eodes  are  used  to  expand  the  researeh  through  simulation  and  analysis 
of  results  obtained  with  ehip  waveform  shaping  and  varying  eode  lengths. 

3.6.2  Workload. 

Many  workload  parameters  exist  and  potentially  impaet  reeeiver  BER.  These 
parameters  inelude; 

•  Signal  strueture  (i.e.,  transmitted  symbol  shape) 

•  Signal  power 

•  Thermal  ehannel  noise 

•  Time  delay  of  interferers  relative  to  the  desired  user 

•  Number  of  ehip  intervals  {Tc)  per  eommunieation  symbol  interval  (T,) 

•  Number  of  eode  periods  per  symbol  interval 

•  Number  of  system  users 

•  Shape  of  ehip  waveforms 

The  transmitted  signals  for  this  researeh  are  baseband  waveforms  eonsisting  of  uniformly 
distributed  random  strings  of  positive  and  negative  ones  that  represent  random  binary 
waveforms.  For  a  given  simulation,  all  users  have  the  same  signal  power.  Thermal 
ehannel  noise  is  modeled  as  AWGN.  Relative  time  delays  for  interfering  users  impaet 
BER  by  ehanging  the  eross-eorrelation  levels  between  the  interfering  eodes  and  the 
desired  user's  eode.  The  number  of  ehips  per  symbol  interval  is  eommonly  related  to 
proeessing  gain,  G,  and  BER  improves  (deereases)  as  G  inereases  [Peterson,  Ziemer,  and 
Borth,  1995].  For  all  simulations  exaetly  one  eode  period  is  eontained  within  eaeh 
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symbol  interval.  In  DS/SSMA  systems  MAI  levels  inerease  as  the  number  of  users 
increases  and  BER  increases.  Finally,  chip  waveform  shape  impacts  code  cross¬ 
correlation  values  such  that  the  dispreading  mixer  output  increases  (or  decreases)  and 
results  in  a  corresponding  increase  (or  decrease)  in  BER. 

5. 7  Factors 

Of  the  parameters  listed  above,  only  the  following  factors  are  varied  for  this 
research;  thermal  noise  power,  time  delay  of  interferers  relative  to  the  desired  user, 
number  of  users,  spreading  code  length,  and  chip  waveform  shape.  Signal  power  is 
accounted  for  in  Eh/No  with  desired  values  achieved  by  holding  Eb  constant  and  varying 
the  noise  power  spectral  density  No.  Eh/No  values  of  0  to  10  dB  are  simulated  for  all 
cases  considered  unless  otherwise  specified.  The  time  delay  of  interferers  relative  to  the 
desired  user  is  varied  between  two  states,  1)  time  delay  equals  zero  for  all  users  in  a 
synchronous  system,  or  2)  time  delays  are  uniformly  distributed  random  variables  in  the 
range  of  [0,  for  an  asynchronous  system.  The  number  of  system  users  directly 
impacts  MAI  levels  and  1)  varies  based  on  availability  of  analytic  approximations  for 
comparison,  or  2)  is  held  constant  at  A  =  9  to  permit  comparison  of  simulated  results  with 
those  presented  by  Cho  and  Eehnert  [Cho  and  Eehnert,  1999].  Spreading  code  length,  N, 
varies  from  31  to  51 1  to  illustrate  the  impact  of  code  length  variation  on  BER.  Finally, 
chip  waveform  shape  varies  based  on  the  analytical  approximations  presented  by  Eehnert, 
and  Kok  and  Do.  Chip  waveform  shapes  considered  include,  rectangular,  Blackman,  and 
Eanczos  as  defined  in  (2.22)  through  (2.25)  [Eehnert,  2002;  Kok  and  Do,  1997].  Figure 
3.2  illustrates  the  normalized  Blackman  and  Eanczos  pulses. 
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Figure  3.2.  Normalized  Chip  Waveform  Shapes,  a)  Blackman  and  b)  Lanczos 

[Kok  and  Do,  1997]. 


3.8  Evaluation  Technique 

All  evaluations  are  based  on  comparison  of  BER  results.  The  BERs  used  are 
either  calculated  based  on  the  BPSK  single-user  BER  expression  of  (2.15),  approximated 
based  on  analytical  expressions  using  (2.14)  and  (2.20)  in  (2.16),  observed  from  previous 
research  as  outlined  in  Chapter  2,  or  estimated  through  simulation  from  this  research. 
Simulation  of  the  physical  system  provides  all  experimental  data  used  for  evaluation. 
Simulations  are  conducted  in  two  phases.  Eirst,  the  code  is  validated  against  known 
analytical  equations  for  a  single  user  and  results  of  previous  research  for  the  case  of 
multiple  access.  Additional  simulations  establish  BER  variation  resulting  from  random 
codes  verses  Gold  codes,  multiple  chip  waveform  shapes,  and  spreading  code  length. 

The  impact  of  factors  is  determined  by  comparing  BER  curves  for  systems  employing 
different  factors  in  the  setup  or  workload.  Using  two  simulation  phases  in  this  research 


3-8 


enables  eode  verifieation  before  exploring  the  impaet  of  additional  faetors  on  BER.  This 
verifieation  supports  the  validity  of  the  new  simulated  eases. 

3.9  Workload 

The  workload  varies  appropriately  for  the  setup  of  eaeh  simulation  phase. 
Workloads  in  all  phases  inelude  baekground  thermal  noise.  Signal  power  levels  used 
provide  adequate  resolution  to  ereate  BER  plots  for  eomparison  with  approximations, 
previous  researeh,  and  new  simulation  results.  Seleetion  of  time  delay  and  proeessing 
gain  refleet  previous  researeh  into  BERs  for  DS/SSMA  systems.  Seleetion  of  three  ehip 
waveform  shapes  enables  quality  eomparison  with  previous  researeh  while  limiting  the 
eomputational  eomplexity. 

3.10  Experimental  Design 

3.10.1  System  Design. 

Matlab®  simulations  are  used  to  produee  BER  vs.  Ei/No  results  for  both 
synchronous  and  asynchronous  systems  with  multiple  users,  using  random  and  Gold 
eodes,  and  implementing  ehip  waveform  shaping.  Eaeh  user  transmits  eommunieation 
waveforms  having  BPSK  baseband  data  modulation  and  eoded  spreading  modulation 
using  the  desired  pulse  shape  over  the  ehip  intervals.  All  signals  are  transmitted 
(reeeived)  with  equal  power  and  noise  power  is  varied  to  aehieve  desired  Ei/No  values. 
Eaeh  ehip  interval  has  duration  Tc  and  eodes  are  periodie  with  length  N  sueh  that  symbol 
duration  Ts  equals  NTc.  All  spreading  pulse  shapes  are  sealed  to  normalize  transmitted 
power  to  that  of  a  system  using  reetangular  ehip  waveforms  (reetangular  waveform 
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shaping  is  the  baseline  for  eomparing  other  waveform  shapes).  Random  spreading  eodes 
are  generated  as  aperiodie,  uniformly  distributed  sequenees  of  positive  and  negative  ones. 
Gold  eodes  are  generated  as  outlined  in  [Peterson,  Ziemer,  and  Borth,  1995]  using  an 
initial  register  state  of  all  ones;  generator  polynomials  [45]8  and  [VSjg  are  used  to  generate 
=  3 1  -length  eodes  and  [  1 02 1  ]  g  and  [  1 46 1  ]  g  are  used  to  generate  =  5 1 1  -length  eodes . 

3.10.2  Code  Validation. 

Code  validation  is  eondueted  in  two  stages.  First,  using  only  reetangular  ehip 
waveforms,  simulation  results  for  1  transmitter  and  A^=  31  length  eodes  (proeessing 
gain)  are  eompared  to  the  single  user  baseline  ealeulated  using  (2.15)  to  validate  eode 
performanee  for  the  simplest  ease.  This  eomparison  is  repeated  using  both  random  and 
Gold  eodes  to  validate  both  eode  implementations  in  the  simulation  eode.  Seeond,  using 
only  random  eodes  and  reetangular  ehip  waveforms,  the  trends  of  simulations  using  K=3 
and  A^  =  3 1  are  eompared  to  findings  of  Gerantiotis  and  Ghaffari  relative  to  the  analytieal 
approximations  for  signal-to-interferenee  plus  noise  (SINR)  in  (2.15)  and  (2.19). 

3.10.3  Waveform  Shaping  and  Code  Length  Tests. 

Following  validation,  simulated  BER  results  are  eompared  for  ehip  waveforms 
having  reetangular,  Blaekman,  and  Lanezos  shapes  as  defined  in  (2.21)  through  (2.25). 
Using  the  verified  simulation  eode,  BER  is  estimated  for  spreading  eodes  eomprised  of 
random  eodes  and  one  full-period  of  the  31 -length  Gold  eodes.  A  eolleetion  of  “best- 
ease”  and  “worst-ease”  Gold  eodes  are  used.  Here,  “best-ease”  Gold  eodes  are  seleeted 
as  the  eombination  of  Gold  eodes  having  the  lowest  eross  eorrelation  values  with  the 
desired  signal.  Eikewise,  “worst-ease”  Gold  eodes  are  seleeted  as  the  eombination  of 
Gold  eodes  having  the  highest  possible  eross  eorrelation  values  with  the  desired  signal. 
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Finally,  to  illustrate  the  impact  of  code  length  on  MAI,  BER  is  simulated  for  systems 
using  A  =  51 1-length  Gold  codes  and  compared  to  aperiodic  random  spreading  codes 
having  A  =  5 1 1  chips  per  symbol  interval. 

3.10.4  Simulation  Length. 

Each  bit  of  data  that  is  coded,  transmitted,  despread,  and  estimated  in  the  receiver 
represents  an  independent  test  of  the  system's  ability  to  accurately  demodulate  the  desired 
user's  data.  Subsequent  bit  transmissions  and  estimations  represent  experimental 
repetitions.  Each  simulation  runs  until  300  bits  are  estimated  in  error.  The  total  number 
of  bits  transmitted,  n,  represents  the  number  of  experimental  trials  in  the  simulation. 
Based  on  BER  values  estimated  in  this  research,  using  a  95%  confidence  interval  and 
running  simulations  until  300  bit  errors  are  detected,  the  sample  mean  observed  through 
simulation  will  vary  from  the  actual  population  mean  due  to  variation  by  an  amount  r  that 
is  at  most  approximately  ±  1 1%  of  the  actual  mean  of  the  simulation  [Canadeo,  2003]. 
Variation  r  is  defined  as 


where  z  =  1 .96  for  the  95%  confidence  interval  and  Pb  is  the  effective  BER  defined  as  the 
ratio  of  the  number  of  bits  in  error  to  n  [Jain,  1991]. 

The  above  analysis  is  based  on  assuming  errors  are  independent  and  normally 
distributed.  The  properties  of  an  AWGN  channel  and  the  central  limit  theorem  indicate 
that  these  errors  should  be  independent  and  normally  distributed.  To  test  these 
assumptions,  control  tests  are  run  with  the  number  of  users  set  at  A  =  5  and  A  =  10  using 
Eh/No  values  of  0  and  5  dB.  Each  test  is  repeated  fifteen  times,  providing  fifteen 
independent  BER  values  for  each  number  of  users  considered  at  a  given  energy  profile. 
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The  experimental  error,  ej,  is  ealeulated  as  the  differenee  between  measured  BER  and  the 
trial  average,  and  this  error  is  plotted  against  the  normal  quantiles,  defined  as 


X,  =4.91 


0.14 


(3.4) 


where 


/-0.5 

(li= - 


(3.5) 


n 

i  is  the  number  of  the  trial  being  plotted,  and  n  is  the  total  number  of  trials.  The  linearity 
of  this  plot  verifies  the  assumption  of  normally  distributed  errors.  Additionally,  the  error, 
ej,  is  plotted  against  the  average  BER  results.  The  lack  of  a  trend  in  the  errors  verifies  the 
independence  of  errors  [Jain,  1991]. 


3.11  Analyze  and  Interpret  Results 

Simulation  results  are  used  to  compile  BER  vs.  Eh/No  curves  for  both  phases  of 
the  experiment.  Relative  benefits  of  tested  factors  are  not  easily  determined  by  direct 
comparison  of  BER  plots  alone  due  to  variation  caused  by  randomness  in  some  of  the 
factors,  e.g.  the  particular  noise  realization.  To  make  statistically  significant 
comparisons,  error  bars  are  included  in  the  BER  curves.  Error  bar  plots  based  on  the 
variation  values  calculated  by  (3.3),  where  the  error  bars  are  plotted  an  amount  r  above 
and  below  the  observed  BER,  provide  bounds  for  actual  BER  mean  values  to  a  95% 
confidence  level  and  allow  more  meaningful  results  through  direct  comparison  [Jain, 
1991].  Eor  BER  values  to  be  “different”  in  a  statistically  significant  manner,  the  error 
bars  of  two  BER  cases  being  compared  cannot  overlap.  Taking  this  into  consideration. 
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direct  comparison  of  BER  curves  yields  meaningful  information  about  the  relative 
performance  enhancements  of  factors  under  consideration. 

Both  phases  of  testing  utilize  direct  comparison  with  error  bars  present  to 
determine  if  results  are  either  significantly  similar  for  validation  or  significantly  different 
for  identifying  performance  enhancements. 

3.12  Summary 

This  chapter  outlines  the  experimental  setup  required  for  1)  validating  the 
DS/SSMA  receiver  simulation  and  2)  comparing  the  relative  performance  enhancements 
of  selected  factors.  The  system's  service  is  the  transmission  of  data  bits  and  metrics  are 
the  BER  vs.  Eh/No  curves  including  error  bars. 

All  system  components  are  simulated,  with  simulations  conducted  in  two  phases. 
The  first  phase  validates  the  analytical  models  of  the  DS/SSMA  receiver  and  the  second 
phase  illustrates  the  relative  impact  of  chip  waveform  shaping,  code  selection,  and  code 
length  on  BER. 

The  analysis  of  results  compares  the  BER  curves  to  expected  theoretical  BER, 
previously  reported  BER  results,  or  other  experimental  results  of  this  research.  Relative 
impact  on  BER  is  the  expected  result  of  these  comparisons. 
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IV.  Results  and  Analysis 


4.1  Validation  Testing 

Simulation  code  is  validated  in  two  stages.  First,  simulated  bit  error  rate  (BER) 
for  a  single  user  in  the  presence  of  noise  is  compared  to  theoretical  BER  calculations 
from  (2.14).  Second,  the  trend  of  simulated  BER  for  multiple  users  is  compared  to 
results  of  Geraniotis  and  Ghaffari  [Geraniotis  and  Ghaffari,  1991]. 


Eigure  4. 1 .  Validation  of  Single  Elser  with  Random  Code  Sequence 
4.1.1  Single  User  Validation. 

Using  only  rectangular  chip  waveforms,  BER  for  a  single  user  with  random 
spreading  codes  is  simulated  for  an  additive  white  Gaussian  noise  (AWGN)  chaimel. 
Error  bars  are  generated  as  outlined  in  Chapter  3  and  reflect  the  95%  confidence  interval 
for  simulated  results.  Simulation  results  for  this  single  user  case  are  shown  in  Eig.  4.1 
along  with  theoretical  BER  calculated  per  (2.14).  As  shown,  the  theoretical  BER  curve 
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lies  within  simulated  error  bar  bounds,  indieating  the  simulated  results  are  eonsistent  with 
theoretical  performance.  Figure  4.2  shows  the  same  relationship  for  the  case  where  Gold 
code  spreading  sequences  are  used  in  place  of  the  random  sequences.  Once  again,  the 
theoretical  result  falls  within  simulated  error  bar  bounds  indicating  simulated  BER  with 
Gold  coded  sequences  is  consistent  with  theoretical  performance. 


Figure  4.2.  Validation  of  Single  User  with  Gold  Code  Sequences 

4.1.2  Multiple  User  Validation. 

Following  the  system  setup  of  Geraniotis  and  Ghaffari,  the  trend  of  simulated 
BER  results  are  compared  to  the  data  from  Table  2.1  [Geraniotis  and  Ghaffari,  1991]. 
The  system  includes  if  =  3  users,  a  spreading  code  of  length  V=  31  chips,  rectangular 
shaped  chip  waveforms,  and  aperiodic  random  spreading  codes.  Simulations  are  run  for 
both  synchronous  and  asynchronous  networks  and  results  are  compared  with  analytical 
approximations  of  (2.19)  and  (2.15),  respectively. 
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Figure  4.3  shows  simulated  results  for  synchronous  users  using  Eb/No  values  of  8, 
10,  12,  14,  and  16  (dB).  Analytic  approximations  from  (2.19)  are  also  shown  for  the 
same  energy  profiles.  Simulation  results  indicate  performance  which  is  poorer  (higher 
BER)  than  predicted  by  the  analytical  approximations.  The  results  presented  here  are 
consistent  with  the  findings  of  Geraniotis  and  Ghaffari  [Geraniotis  and  Ghaffari,  1991]. 
Given  the  analytic  approximation  falls  well  outside  the  error  bar  bounds  for  all  values  of 
Et/No  considered,  these  results  are  statistically  significant  at  the  95%  confidence  level. 


Figure  4.3.  Simulated  vs.  Analytic  BER  for  Synchronous  DS/SSMA  Systems  Using 

Random  Code  Sequences. 

A  comparison  of  calculated  BER  from  (2.15)  with  simulation  results  of  an 
asynchronous  system,  using  Eft/Ao  values  of  8,  10,  12,  14,  and  16  (dB),  is  shown  in 
Fig.  4.4.  Once  again,  the  data  indicates  that  simulated  BER  under  performs  (higher  BER) 
the  approximations  and  are  consistent  with  Geraniotis  and  Ghaffari's  findings  that  the 
analytical  approximation  is  “optimistic”  for  the  asynchronous  case  [Geraniotis  and 
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Ghaffari,  1991].  Again,  the  error  bars  add  statistical  significance  to  this  result  at  the  95% 
confidence  level. 


Figure  4.4.  Simulated  vs.  Analytical  BER  ^or  Asynchronous  DS/SSMA  Systems  Using 

Random  Code  Sequences. 

4.1.3  Error  Analysis 

Error  analysis  was  conducted  to  verify  the  accuracy  of  the  synchronous 
simulation  model.  Using  only  rectangular  pulse  shapes  and  random  code  sequences, 
simulations  for  A  =  5  and  10  users  at  Eh/No  values  of  0  and  5  (dB)  are  repeated  fifteen 
times.  Experimental  error,  cj,  is  calculated  as  the  difference  between  the  observed  BER 
for  each  trial  and  the  average  of  fifteen  trials.  Values  of  cj  are  plotted  against  the  normal 
quantiles  as  determined  by  (3.4)  and  (3.5)  in  fig.  4.5,  fig.  4.6,  fig.  4.7,  and  fig.  4.8.  The 
nearly  linear  nature  of  these  plots  verifies  the  errors  are  normally  distributed  [Jain,  1991]. 
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Figure  4.5.  Quantile-Quantile  Plot  for  =  5  Users  and  Ei/No  =  0  (dB). 


Figure  4.6.  Quantile-Quantile  Plot  for  if  =  5  Users  and  Et/No  =  5  (dB). 
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Figure  4.7.  Quantile-Quantile  Plot  for  =  10  Users  and  Ei/No  =  0  (dB). 
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Figure  4.8.  Quantile-Quantile  Plot  for  =  10  Users  and  Eb/No  =  5  (dB). 
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Figure  4.9  and  Fig.  4.10  show  the  residual  quantity,  ej,  plotted  about  zero,  whieh 
represents  the  average  BER  from  the  trials.  All  residual  errors  are  uniformly  distributed 
about  the  trial  mean  indieating  a  lack  of  bias.  Note;  The  first  fifteen  points  are  for 
Eh/No  =  0,  and  the  next  fifteen  points  are  for  Et/No  =  5.0  in  Fig.  4.9  and  Fig.  4.10.  As 
Eh/No  increases  BER  decreases  and  residual  variation  from  the  mean  also  decreases.  This 
is  consistent  with  the  prediction  that  the  errors  would  fall  within  ±  1 1%  of  the  mean  BER 
[Jain,  1991]. 


0.015 


0.01 


0.005 


Residual 

Value 


-0.005 


-0.01 


-0.015 


Residual 

—  Average  BER 


♦  ♦  ♦  ♦ 

♦  ♦  ♦ 


10  15  20 

Trial  Number 


25 


30 


Eigure  4.9.  Residual  Plot  for  A  =  5  Users. 
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Figure  4.10.  Residual  Plot  for  Ai  =  10  Users. 

4.1.4  Conclusions. 

The  two  stage  validation  tests  illustrate  that  the  fundamental  system  model  and 
eoding  in  the  Matlab®  simulation  are  sufficiently  accurate  to  perform  comparative 
analyses.  Error  analysis  plots  indicate  the  system  model  as  simulated  is  free  from 
undesired  bias.  The  results  of  these  validation  and  error  analysis  tests  provide  a  solid 
basis  for  establishing  model  accuracy  and  support  the  experimental  results  of  this 
research. 


4.2  Impact  of  Gold  Codes 

Previous  research  into  the  impact  of  multiple  access  interference  (MAI)  on 
DS/SSMA  system  BER  has  focused  primarily  on  systems  using  aperiodic  random 
spreading  code  sequences,  or  m-Sequences  [Cho  and  Eehnert,  1999;  Geraniotis  and 
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Ghaffari,  1991].  Results  of  simulation  for  a  DS/SSMA  system  with  A  =  9  users,  A  =  31 
length  eodes,  using  reetangular  ehip  waveforms  and  employing  Gold  eode  sequenees  is 
eompared  to  BER  results  for  aperiodie  random  spreading  eodes  for  both  synchronous  and 
asynchronous  systems. 

Figure  4.1 1  shows  BER  eurves  for  both  best-ease  and  worst-ease  Gold  eodes 
(relative  to  random  eodes)  for  the  synchronous  system.  The  best-ease  Gold  eoded  system 
far  outperforms  the  randomly  eoded  system,  while  the  worst-ease  Gold  eoded  system 
under  performs  the  randomly  eoded  system.  Error  bars  are  ineluded  in  Fig.  4.1 1  and 
indieate  the  differenees  noted  are  statistieally  signifieant  for  all  Et/No  values  tested. 

Figure  4.12  eompares  BER  results  for  best-ease  and  worst-ease  Gold  eodes 
(relative  to  random  eodes)  for  the  asynchronous  system.  Again,  the  best-ease  Gold  eoded 
system  outperforms  the  randomly  eoded  system  and  the  worst-ease  Gold  eoded  system 
under  performs  the  randomly  eoded  system.  The  error  bars  indieate  this  performanee 
improvement  is  signifieant  for  Eh/No  values  greater  than  1 .0  (dB).  However,  the 
performanee  differenees  due  to  Gold  eoding  noted  for  the  asynchronous  system  are  less 
than  observed  for  the  synchronous  system.  These  results  illustrate  that  Gold  eode 
oharaeteristies,  i.e.,  the  eross-eorrelation  and  autoeorrelation  values  for  a  given  family, 
are  well  suited  for  the  synchronous  ease  eonsidered  [Peterson,  Ziemer,  and  Borth,  1995]. 
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Figure  4.1 1.  Synchronous  Network;  Gold  vs.  Randomly  Coded  System 


Figure  4.12.  Asynchronous  Network:  Gold  vs.  Randomly  Coded  System 
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4. 3  Pulse  Shaping  Results 

The  effects  of  chip  waveform  shaping  on  synchronous  and  asynchronous 
DS/SSMA  system  performance  are  characterized  via  simulation.  All  simulations 
contained  A  =  9  users  with  A  =  31 -length  spreading  codes  to  permit  comparison  of  results 
obtained  by  Cho  and  Lehnert  [Cho  and  Lehnert,  1999].  BER  results  are  obtained  using 
randomly  coded  sequences,  best-case  and  worst-case  Gold  coded  sequences  and 
rectangular,  Blackman,  and  Lanczos  shaped  chip  waveforms. 

4.3.1  Synchronous  System  Performance 

As  shown  in  Fig.  4.13,  Fig.  4.14,  and  Fig.  4.15,  simulation  results  for 
synchronous  DS/SSMA  systems  show  very  little  improvement  in  BER  performance  as 
the  chip  waveform  shape  changes;  the  resultant  error  bars  overlap  which  indicates  the 
differences  shown  are  not  statistically  significant. 


4-11 


Figure  4.13.  BER  for  Synchronous  Randomly  Coded  DS/SSMA  System  with  K=9 
Users,  A  =  3 1  Length  Codes  and  Chip  Waveform  Shaping. 


Figure  4.14.  BER  for  Synchronous  Best-Case  Gold  Coded  DS/SSMA  System  with 
K  =  9  Users,  N  =  31 -Length  Codes  and  Chip  Waveform  Shaping. 
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Figure  4.15.  BER  for  Synchronous  Worst-Case  Gold  Coded  DS/SSMA  System  with 
K=9  Users,  A  =  3 1  Length  Codes  and  Chip  Waveform  Shaping. 

4.3.2  Asynchronous  System  Performance. 

Results  for  the  asynchronous  system  using  randomly  coded  sequences  are  shown 
in  Fig.  4.16.  Clearly,  the  BER  resulting  from  the  Lanczos  waveform  shape  is  improved 
relative  to  the  rectangular  waveform  shape;  an  improvement  of  approximately  6.0  dB  is 
indicated.  Results  presented  in  Fig.  4.17  for  the  best-case  Gold  coded  asynchronous 
system  also  indicate  improvement  when  the  Lanczos  waveform  shape  is  introduced;  an 
improvement  of  approximately  4.5  dB  is  indicated.  As  seen  in  Fig.  4.18,  when  the  worst- 
case  Gold  coded  sequences  are  used  in  an  asynchronous  system,  an  improvement  for 
Lanczos  waveform  shapes  of  approximately  7.5  (dB)  is  realized. 
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Figure  4.16.  BER  ^or  Asynchronous  Randomly  Coded  DS/SSMA  System  with  A  =  9 
Users,  A  =  3 1  Length  Codes  and  Chip  Waveform  Shaping. 


Figure  4.17.  BER  for  Asynchronous  Best-Case  Gold  Coded  DS/SSMA  System  with 
K=9  Users,  A  =  3 1  Length  Codes  and  Chip  Waveform  Shaping. 
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Figure  4.18.  BER  ^ox  Asynchronous  Worst-Case  Gold  Coded  DS/SSMA  System  with 
K=9  Users,  A  =  3 1  Length  Codes  and  Chip  Waveform  Shaping. 


4.4  Cross-Correlation  Analysis 

The  particular  fixed  propagation  delays  associated  with  asynchronous  results 
presented  in  Section  4.3  are  201,  263,  273,  279,  309,  310,  171,  and  339  samples  for  users 
2  through  9,  respectively.  These  delays  represent  Delay  Profile  #1.  Cross-correlation 
values  between  delayed  interferers  and  the  desired  signal  are  calculated  for  each  chip 
waveform  shape  considered  and  are  listed  in  Table  4.1  along  with  the  sum  of  the 
interferer  cross-correlations.  The  sum  of  cross-correlation  values  for  the  best-case  Gold 
coded  case  progressively  decreases  when  going  from  rectangular,  to  Blackman,  to 
Lanczos  pulse  shaping.  Cumulative  cross-correlation  values  obtained  from  Lanczos 
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pulse  shaping  is  a  factor  of  6.25  lower  than  that  of  the  rectangular  pulse  shape.  This  can 
be  mapped  to  a  corresponding  performance  gain  of  7.96  dB  as  obtained  from: 


Processing  Gain  =  10  x  logjg 


Correlation  of  Interest^ 
Reference  Correlation 


(4.1) 


This  calculated  processing  gain  is  greater  than  that  reflected  in  BER  curves;  Lanczos 
waveform  shaping  provides  an  approximate  4.5  dB  improvement  over  rectangular 
waveform  shaping.  Systems  using  worst-case  Gold  coded  sequences  show  a  factor  of  2.1 
decrease  in  cumulative  cross-correlation  values  between  rectangular  and  Lanczos 
waveform  shapes,  corresponding  to  an  approximate  3.2  dB  improvement;  BER  curves  of 
Eig.  4.18  indicate  much  greater  improvement. 


Table  4. 1 .  Cross-Correlation  Between  Desired  User  and  Interfering  Users  Employing 
Best-Case  Gold  Codes  and  Chip  Waveform  Shaping  for  Delay  Profile  #1. 


User  # 

2 

3 

4 

5 

6 

7 

8 

9 

Delay 

201 

263 

273 

279 

309 

310 

171 

339 

Sum 

Cross-Correlation  Value  Multiplied  by  341 

Rectangular 

-11 

-11 

61 

-11 

77 

5 

29 

5 

144 

Blackman 

-4.83 

-10.13 

54.82 

-2.12 

70.91 

-7.89 

-1.78 

-7.89 

91.09 

Lanczos 

-4.69 

-2.27 

11.04 

-7.20 

15.86 

3.38 

3.54 

3.38 

23.05 

Each  multiple  access  interferer  is  modeled  as  having  Ns  x  Nc  (number  of  samples 
per  symbol  times  the  number  of  chips  per  symbol)  possible  delay  values;  for  parameters 
used  here  there  are  341  possible  delay  values  and  {K-\)  =  8  total  multiple  access 
interferers  in  the  system.  Therefore,  there  are  341  =  1.8  x  10  possible  combinations  of 
delay  values  for  the  system  under  consideration  [Leon-Garcia,  1994].  To  account  for  the 
large  number  of  possible  delay  profiles  without  conducting  exhaustive  Monte  Carlo 
simulations,  the  average  cross-correlation  value  for  each  interfering  code  using 
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rectangular  waveform  shapes,  and  the  cumulative  of  those  averages,  as  listed  in 
Table  4.2,  are  considered. 


Table  4.2  Average  Cross-Correlation  Between  Desired  User  and  Interfering  Users 
Employing  Best-Case  Gold  Codes  for  All  Possible  Delays. 


User  # 

2 

3 

4 

5 

6 

7 

8 

9 

Total 

-1.243 

0.1775 

0.1775 

-1.243 

-1.243 

0.1775 

-1.2453 

0.1775 

-4.2643 

For  simulation  purposes,  specific  fixed  propagation  delay  values  for  the 
interferers  were  selected  such  that  the  cumulative  cross-correlation  values  obtained  in  the 


simulation  closely  matched  the  magnitude  of  the  values  listed  in  Table  4.2.  The  resulting 
delays  for  users  2  through  9  are  listed  in  Table  4.3  (Delay  Profile  #2)  and  produce  a 
cumulative  cross-correlation  value  of  8/341  for  rectangular  chip  waveforms.  The  BER 
performance  is  again  simulated  using  Delay  Profile  #2  with  results  shown  in  Fig.  4.19, 
Fig.  4.20,  and  Fig.  4.21  for  random,  best-case  Gold  and  worst-case  Gold  coded 
sequences,  respectively. 


Table  4.3  Delay  Profile  #2  Based  on  Achieving  Near  Average  Cross-Correlation 
Between  Desired  User  and  Interfering  Users  Employing  Best-Case  Gold  Codes  and 

Chip  Waveform  Shaping. 


User  # 

2 

3 

4 

5 

6 

7 

8 

9 

Delay 

17 

134 

266 

93 

329 

207 

155 

13 

Sum 

Cross-Correlation  Value  Multiplied  by  341 

Rectangular 

-3 

5 

5 

-3 

-3 

5 

-3 

5 

8 

Blackman 

6.11 

-7.89 

-7.89 

6.11 

-9.94 

-7.89 

-9.94 

-7.89 

-39.25 

Fanczos 

3.62 

3.38 

3.38 

3.62 

-1.97 

3.38 

-1.97 

3.38 

16.83 
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Figure  4.19.  BER  for  Randomly  Coded  DS/SSMA  System  with  Seleeted  Delays; 
K=9  Users  and  A=  3 1  Length  Codes. 


Figure  4.20.  BER  for  Best-Case  Gold  Coded  DS/SSMA  System  with  Selected  Delays; 

K=9  Users  and  A=  3 1  Length  Codes. 
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Figure  4.21.  BER  for  Worst-Case  Gold  Cot/e(i  DS/SSMA  System  with  Seleeted 
Delays;  K=9  Users  and  A  =  3 1  Length  Codes. 

Using  the  fixed  Delay  Profile  #2  deseribed  by  values  in  Table  4.3,  the  BER  of  a 
best-ease  Gold  eoded  system  shows  BER  improvement  of  0.0  and  3.5  dB  for  Blackman 
and  Lanczos  pulse  shaping,  respectively.  These  improvements  are  not  consistent  with  the 
processing  gain  of  -6.91  and  -3.22  dB  obtained  by  using  the  cumulative  average  cross¬ 
correlation  values  in  Table  4.3,  indicating  an  error  in  either  the  prediction  or  the 
simulation. 

To  isolate  the  potential  source  of  error,  a  more  detailed  analysis  of  the  test  statistic 
involved  in  this  process  is  conducted.  The  total  test  statistic  z  (the  correlator  output  for 
signal,  MAI  and  noise  present)  is  defined  as; 

9 

^  =  (4-2) 

j=2 
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where  z;^  is  the  autoeorrelation  response  between  user  one's  transmitted  (spread)  symbol 
with  user  one's  spreading  eode  (the  desired  user),  the  zij  terms  in  the  summation  are  the 
eross-eorrelation  response  between  user  one's  spreading  eode  and  the user's  transmitted 
symbol,  and  z„  is  the  noise  component  out  of  the  correlator.  Table  4.4  lists  the  mean  and 
variance  of  individual  test  statistics  for  an  Et/No  value  of  3.0  dB  using  Delay  Profde  #2, 
where  zn  is  the  summation  term  from  (4.2)  representing  the  summation  of  cross¬ 
correlations  of  each  interfering  user  with  user  one.  The  zn  term  is  the  determining  factor 
in  predicting  MAI  effects.  The  variances  of  the  test  statistic  components  indicate  the 
amount  of  power  associated  with  each  term. 

Table  4.4.  Test  Statistic  Mean  and  Variance  Values 


Zll 

Zn 

'^n 

M 

(7^(X10^) 

f7^(xl0^) 

fi 

(xlO^) 

Rectangular 

6.019 

1.163 

-0.177 

6.325 

0.005 

2.880 

Blackman 

0.597 

1.163 

-0.300 

6.227 

1.651 

2.890 

Eanczos 

0.424 

1.163 

-0.036 

0.292 

-0.098 

2.898 

For  the  case  where  Blackman  shaped  chip  waveforms  are  used,  the  variance  in  the 
MAI  term  (z;,)  is  not  significantly  different  from  the  rectangular  shaped  chip  waveforms. 
This  trend  is  consistent  with  BER  data  plotted  in  Fig.  4.20  which  shows  little  to  no 
variation  between  Blackman  and  rectangular  pulse  shaping.  As  shown  in  Table  4.4,  the 
variance  in  the  MAI  term  for  the  Lanczos  shaped  chip  waveforms  is  approximately  one 
order-of-magnitude  lower  than  both  the  rectangular  and  Blackman  case,  indicating 
significantly  better  BER  performance  should  be  realized  for  Eanczos  shaped  waveforms 
relative  to  either  rectangular  or  Blackman  shaped  chip  waveforms.  This  trend  is  also 
consistent  with  the  BER  results  in  Fig.  4.20. 


4-20 


Due  to  the  temporal  overlap  of  the  desired  user's  symbol  interval  with  two 
portions  of  eaeh  interfering  users'  symbol  intervals,  there  exist  four  possible  data  value 
combinations  for  each  asynchronous  interfering  user  in  the  despreading  process;  the 
interferer  can  have  a  value  of  either  one  or  zero  for  both  symbols,  a  one  then  a  zero,  or  a 
zero  then  a  one.  These  data  combinations  effect  the  cross-correlation  between  the  desired 
user's  spreading  code  and  the  delayed  interfering  user's  spreading  code.  Table  4.5  and 
Table  4.6  list  the  mean  and  variance  of  the  cross-correlations  for  users  2  through  9  when 
these  data  effects  are  included.  The  summation  of  the  variances  of  the  cross-correlations 
with  data  effects  results  in  6.330  x  10^,  6.195  x  10^,  and  0.294  x  10^  for  rectangular, 
Blackman,  and  Lanczos  chip  waveforms,  respectively.  These  values  indicate  that  for  the 
specific  delays  listed  in  Table  4.3  the  systems  with  rectangular  and  Blackman  chip 
waveforms  should  perform  similarly,  while  the  system  with  Lanczos  chip  waveforms 
should  significantly  outperform  the  other  two.  The  simulation  results  in  Fig.  4.20  reflect 
this  prediction. 

Table  4.5.  Mean  and  Variance  of  Cross-Correlation  with  Data  Effects  for 


Users  2  through  5  using  Delay  Profile  #2 


User  # 

2 

3 

4 

5 

Delay 

17 

134 

266 

93 

xlO^ 

n 

xlO^ 

n 

xlO^ 

£7^  XIO^ 

Rectangular 

0.019 

0.029 

0.053 

0.192 

0.074 

3.815 

-0.026 

0.009 

Blackman 

0.222 

0.035 

-0.050 

0.800 

0.063 

2.502 

-0.042 

0.022 

Lanczos 

-0.024 

0.019 

0.042 

0.006 

-0.008 

0.157 

0.009 

0.009 
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Table  4.6.  Mean  and  Variance  of  Cross-Correlation  with  Data  Effects  for 
Users  6  through  9  using  Delay  Profile  #2 


User  # 

6 

7 

8 

9 

Delay 

329 

207 

155 

13 

xlO^ 

CT  xlO^ 

0-^  xlO^ 

xlO^ 

Reetangular 

0.011 

0.362 

0.029 

1.638 

-0.013 

0.186 

-0.026 

0.098 

Blaekman 

0.016 

0.515 

-0.034 

1.505 

0.177 

0.505 

0.026 

0.310 

Lanezos 

-0.003 

0.024 

0.025 

0.052 

0.006 

0.022 

0.074 

0.006 

The  cross-correlation  statistics  for  the  system  design  resulting  in  Fig.  4.17  BER 
curves  are  listed  in  Table  4.7.  and  Table  4.8.  The  summation  of  cross-correlation 
variances  including  data  effects  results  in  1 1.185  x  10  ,  5.537  x  10  and  0.953  x  10  for 
rectangular,  Blackman,  and  Eanczos  chip  shapes,  respectively.  These  relative  variance 
values  are  consistent  with  the  significant  improvement  displayed  in  Fig.  4.17.  for 
Blackman  and  Eanczos  chip  waveforms  over  rectangular  chip  waveforms  for  Delay 


Profile  #1. 

Table  4.7.  Mean  and  Variance  of  Cross  Correlation  with  Data  Effects 
for  Users  2  through  5  using  Delay  Profile  #1 


User  # 

2 

3 

4 

5 

Delay 

201 

263 

273 

279 

xlO^ 

xlO^ 

0-^  xlO^ 

a  xlO^ 

Reetangular 

0.162 

0.071 

-0.746 

0.441 

-7.018 

2.890 

0.027 

0.536 

Blaekman 

0.306 

0.025 

0.214 

0.535 

-0.495 

1.742 

0.045 

0.082 

Lanezos 

0.000 

0.012 

-0.049 

0.024 

0.238 

0.128 

-0.067 

0.525 

Table  4.8.  Mean  and  Variance  of  Cross  Correlation  with  Data  Effects 
for  Users  2  through  5  using  Delay  Profile  #1 


User  # 

6 

7 

8 

9 

Delay 

309 

310 

171 

339 

tExlO^ 

ff^xlO^ 

tExlO^ 

ff^xlO^ 

Reetangular 

-3.673 

2.985 

0.198 

0.096 

-0.251 

4.111 

0.440 

0.055 

Blaekman 

-1.495 

2.791 

-0.142 

0.296 

-0.081 

0.003 

-0.309 

0.062 

Lanezos 

0.371 

0.128 

0.026 

0.006 

-0.105 

0.112 

0.165 

0.017 

4-22 


4.5  Code  Length  Analysis 

The  BER  performance  of  a  system  using  A^=  31 -length  spreading  codes  is 
compared  to  the  BER  performance  of  a  system  using  A^=  51 1-length  spreading  codes  as 
shown  in  Eig.  4.22.  In  this  case,  data  is  presented  for  synchronous  systems  with  K=9 
users,  random  and  best-case  Gold  coded  sequences,  and  rectangular  waveform  shapes. 

All  cases  use  a  fixed  symbol  duration  and  one  full  period  of  the  spreading  code  per 
symbol  duration.  Data  in  Eig.  4.22  reflects  an  improvement  of  approximately  7.5  dB  and 
0.5  dB  for  random  and  Gold  codes,  respectively,  as  code  length  increases  from  31  to  51 1. 
Eor  the  random  case,  although  appreciable,  this  improvement  is  less  than  the  theoretical 
processing  gain  improvement  of  12.7  dB  expected  by  increasing  N  from  31  to  5 1 1 .  The 
Gold  coded  case  shows  some  improvement,  but  falls  far  short  of  the  expected  12.7  dB 
gain,  and  is  discussed  below. 

BER  results  for  5 1 1 -length  randomly  and  Gold  coded  results  for  =  9  are 
compared  to  the  single  user  communications  performance  calculated  from  (2.15)  as 
shown  in  Eig.  4.23.  Implementing  51 1-length  spreading  codes  increases  the  9  user 
system  performance  (decreases  BER)  to  a  level  that  is  almost  equivalent  to  the  single  user 
case  for  randomly  coded  waveforms.  Eor  Gold  coded  systems,  BER  improves  to  a  level 
that  is  not  statistically  different  than  the  single  user  communications  performance  (no 
MAI  present).  Due  to  the  presence  of  AWGN,  system  BER  performance  is  limited 
(cannot  decrease  below)  the  single-user  communications  performance.  Marginal 
performance  gain  is  observed  when  increasing  the  Gold  code  length  from  3 1  to  5 1 1  given 
performance  cannot  be  any  better  than  the  case  where  there  is  no  MAI  present. 
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Figure  4.22.  BER  for  Randomly  and  Gold  Coded  Synchronous  DS/SSMA  Systems 
with  Fixed  Symbol  Duration;  K=9  Users  and  A  =  3 1  and  511  Length  Codes  using 
One  Full  Period  of  the  Spreading  Code  per  Symbol  Duration. 


Figure  4.23.  Comparison  of  Single  User  Performanee  to  Randomly  and  Gold  Coded 
Synchronous  DS/SSMA  Systems;  K=9  Users  and  A  =  5 1 1  Length  Codes. 
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V.  Conclusions 


5.1  Research  Contributions 

This  research  provides  characterization  of  multiple  access  interference  (MAI)  for 
direct-sequence  spread-spectrum  multiple  access  (DS/SSMA)  systems  through  variation 
of  several  factors,  including  the  use  of  random  versus  Gold  spreading  codes,  chip 
waveform  shape  selection,  and  variation  in  spreading  code  length.  Results  provide 
insight  into  the  applicability  of  previously  published  approximations  for  bit  error  rate 
(BER)  and  the  impact  that  each  system  factor  has  on  overall  performance.  This  research 
also  provides  a  modeling  and  simulation  tool  capable  of  supporting  future  research  into 
the  factors  effecting  DS/SSMA  system  performance.  This  simulated  model  enables  a 
more  realistic  implementation  of  potential  system  configurations  to  verify  future 
analytical  approximations. 

5.2  Summary  of  Findings 

5.2.1  Random  versus  Gold  Spreading  Codes 

For  DS/SSMA  systems  employing  Gold  spreading  codes,  performance  is  shown 
to  be  greater  than  similar  systems  employing  random  spreading  codes.  This 
improvement  is  most  noticeable  at  Et/No  values  greater  than  1 .0  dB,  the  point  at  which 
the  effect  of  MAI  on  BER  is  large  enough  to  dominate  the  effect  of  additive  white 
Gaussian  noise  (AWGN)  on  BER.  Additionally,  the  improvement  due  to  Gold  coding  is 
more  pronounced  in  synchronous  systems,  indicating  that  the  power  of  Gold  coding  is 
based  on  the  synchronous  characteristics  of  the  codes. 
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5.2.2  Chip  Waveform  Shaping 

Improvement  in  BER  due  to  ehip  waveform  shaping  is  shown  to  be  mostly 
limited  to  asynchronous  systems;  changes  in  spreading  code  cross-correlation  statistics 
(relative  to  the  synchronous  case)  occur  due  to  symbol  boundary  misalignment  in  the 
asynchronous  case.  For  the  case  of  fixed,  predetermined  propagation  delay  profiles  for 
the  synchronous  users,  the  BER  improves  (or  degrades)  by  an  amount  that  is  in  part 
dictated  by  the  improvement  (or  degradation)  in  variance  of  symbol  cross-correlation 
which  occurs  when  introducing  an  alternate  chip  waveform  shape;  improved  cross¬ 
correlation  variance  here  implies  lowers  MAI  levels  and  corresponding  lower  BERs  as 
the  chip  waveform  shape  is  changed.  The  demonstrated  improvement  in  this  work  is  an 
effective  processing  gain  of  between  3.5  and  4.5  dB  for  a  system  using  Eanczos  chip 
waveforms,  relative  to  identical  systems  using  rectangular  chip  waveforms  in  an 
asynchronous  network  with  best-case  Gold  codes.  This  improvement  was  found  to  be 
dependent  on  the  variance  of  the  cross-correlation  between  spreading  codes  when 
considering  data  effects. 

5.2.3  Spreading  Code  Length 

An  increase  in  spreading  code  length  from  3 1  chips  to  5 1 1  chips  provides 
significant  improvement  in  BER  for  a  synchronous  DS/SSMA  system  employing  random 
spreading  codes.  For  a  given  BER,  a  reduction  in  required  Et/No  of  approximately  7.5  dB 
is  demonstrated,  falling  below  the  expected  performance  improvement  (reduction)  of 
10xlogio(E  =  511/31)-  12.7  dB.  In  cases  where  total  interfering  power  is  dominated  by 
MAI  and  not  by  AWGN,  increasing  code  length  by  a  factor  off  effectively  reduces  the 
dominant  MAI  power  by  an  amount  proportional  to  ML  and  decreases  (improves)  overall 


5-2 


BER.  For  Gold  coded  systems  the  performanee  improvement  was  only  0.5  dB  for  the 
same  eode  length  inerease.  This  improvement  was  shown  to  deerease  BER  to  the  level  of 
single  user  performanee,  whieh  is  the  theoretical  best  performance  in  the  presence  of 
AWGN. 

5.5  Recommendations  for  Future  Research 

5.3.1  M-Ary  Data  and/or  Spreading  Modulation 

This  researeh  only  simulated  binary  modulation  and  spreading  teehniques,  i.e., 
binary  phase  shift  keyed  (BPSK)  data  and  spreading  modulations.  DS/SSMA  systems 
ean  be  implemented  using  other  M-ary  modulations,  e.g.,  BPSK  data  modulation  ean  be 
eombined  with  quadrature  phase  shift  keyed  (QPSK)  spreading  modulation  as  done  in  the 
eellular  IS-95  eommunication  system.  The  impaet  on  BER  performance  using  alternate 
modulation  teehniques  eould  be  easily  simulated  with  only  minor  modification  of 
simulation  code  generated  under  this  researeh. 

5.3.2  Full  Monte  Carlo  Simulations 

Given  the  randomness  assoeiated  with  several  faetors  in  a  DS/SSMA  system,  i.e., 
random  multiple  aecess  coding,  asynchronous  user  propagation  delays,  etc.,  several 
variables  were  fixed  for  simulations  in  this  work  sueh  that  they  provided  “average” 
operating  charaeteristics  for  the  simulation.  A  full  Monte  Carlo  simulation  of  the  system 
variables  eould  be  used  to  more  fully  eharaeterize  BER  performanee  in  the  presenee  of 
MAI. 
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5.5.5  Near-Far  Comparisons 

Scenarios  could  be  considered  where  interfering  users  are  transmitting  at  locations 
either  “nearer”  to  or  “farther”  from  the  receiver  of  interest.  Simulation  of  these  scenarios 
would  address  the  so  called  “near-far”  problem  which  was  not  addressed  as  part  of  this 
research  (in  this  work  all  interfering  signals  were  assumed  to  be  received  at  the  same 
power  level  as  the  desired  user).  In  near-far  cases,  the  MAI  power  received  varies  from 
user-to-user  relative  to  the  desired  user  and  this  variation  influences  the  impact  of  other 
factors  on  BER. 

5.3.4  Characterization  Using  Delay  Profile  Variation 

For  this  research,  variances  in  cross-correlation  (which  directly  correlate  to  MAI 
power  levels)  were  used  to  induce  the  interaction  of  various  spreading  codes  on  expected 
performance  improvement.  For  the  asynchronous  network  case,  these  MAI  power  levels 
were  obtained  by  using  a  specific  propagation  delay  profile  (collection  of  pre-assigned 
delay  values)  for  the  simulations.  As  revealed  throughout  the  preliminary  stages  of  this 
research,  there  exist  specific  propagation  delay  profiles  whereby  a  change  in  chip 
waveform  shape  on  interfering  users  did  not  improve  BER  performance  relative  to  the 
rectangular  chip  waveforms.  While  average  performance  improvement  provides  valuable 
academic  insight,  accurately  characterizing  performance  for  delay  profile  variation 
(which  can  occur  with  changes  in  system  topology)  may  be  important,  especially  for 
actual  communication  systems. 
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Appendix  A.  Simulation  Code 


Matlab®  simulation  code  is  presented  in  this  appendix.  The  primary  funetion  used  for  31- 
length  eodes  was  ber_shape_async_all31_2.m  and  for  51 1-length  eodes  was 
ber_shape_syne_BR5 1 1  .m. 

A .  1  ber_shape_async_aU3 1_2.  m 

function  [berb,  berr,  berw]  =  ber_shape_async_all31_norm2(eb,numbsig, samps, delay, shape) 

% 

%  [berb  berr  berw]  =  ber_shape_async_rand31_2(eb,numbsig,samps,delay,shape) 

% 

%  Change  the  random  spreading  codes  to  be  aperiodic. 

% 

%  Experimentally  determine  the  BER  for 
%  numbsig  users  at  each  of  the  eb  values 
%  of  Eb/No  (dB).  All  signals  are  asynchronous. 

% 

%  Inputs:  eb  =  Eb/No  (dB) 

%  numbsig  =  number  of  signals  (including  desired  user) 

%  samps  =  samples  per  chip  interval 
%  srand  =  random  codes 

%  delay  =  number  of  smaples  time  delay  for  each  user  (from  0  - 
%  total  samples  per  symbol  interval) 

%  shape  =  desired  shape  for  the  chip  interval 
%  1  =  rectangular 

%  2  =  half  sine 

%  3  =  raised  cosine 

%  4  =  blackman 

%  5  =  raised  cosine  +  1 

% 

%  Returns 

%  berb  =  BER  using  best  case  Gold  codes 
%  berw  =  BER  using  worst  case  Gold  codes 
%  berr  =  BER  using  random  codes  defined  by  srand 

% 

% 

%  Written  by  Matthew  Glen  5  Dec  03 

%%%%  Spread  Spectrum  MAI  characterization  project 
%%%%%  Create  m-seq  needed  for  Gold  codes 

bits  =  1 ; 
slen  =  31; 
nsamps  =100; 
t  =  0:samps-l; 

if  shape  ==  1 
wf=  ones(l, samps); 
elseif  shape  ==  2 

wf  =  sin(pi*t/(length(t)-l)); 
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elseif  shape  ==  3 

wf  =  sqrt(2/3)*(l  -  cos(2*pi*t/(length(t)-l))); 
elseif  shape  ==  4 

wf  =  0.42  -  0.5*eos(2*pi*t/(length(t)-l)); 

elseif  shape  ==  5 
t2  =  t+10; 

wf  =  (siii(2*pi*t2)./(2*pi*t2)).^2; 
end 

%  gpolyl  =  [1  0000  1  000  1];  %%%  1021 

%  gpoly2  =  [1  1  00  1  1  000  1];  %%%  1461 

%  init  =  [111111111];  %%%  initial  register  state 

%  mseq  =  [zeros(2,5 1 1)];  %%%  ereate  mseq  veetor  of  zeros 

%  mseq(l,:)  =  mSeqGen(init, gpolyl);  %%%  generate  the  first  mseq 
%  mseq(2,:)  =  mSeqGen(init,gpoly2);  %%%  generate  the  2nd  mseq 

gpolyl  =[100101];  %%%  45 

gpoly2  =  [111101];  %%%  75 

init=  [11111];  %%%  initial  register  state 

mseq  =  [zeros(2,3 1)];  %%%  ereate  mseq  veetor  of  zeros 

mseq(l,:)  =  mSeqGen(init, gpolyl);  %%%  generate  the  first  mseq 

mseq(2,:)  =  mSeqGen(init,gpoly2);  %%%  generate  the  2nd  mseq 

gold  =  gold_gen(mseq(l,:),  mseq(2,:)).*(-2)  +1;  %%  see  whieh  direetion  they  are  lined  up 
[row  eol]  =  size(gold); 

%%  Run  test  until  endbit  number  of  errors  is  reaehed 
endbit  =  300; 

%  Define  the  best  and  worst  eombinations  of  51 1-length  gold  eodes 
bestorder  =  [1  3:33];  %  best  eodes 

worstorder  =  [2  33  30  5  7  10  17  31  29  27  26  24  22  21  13  4  1];  %  worst  eodes 
%%%  Spreading  eodes 

best(l:numbsig,:)  =  gold(bestorder(l:numbsig),:); 
worst(l:numbsig,:)  =  gold(worstorder(l:numbsig),:); 

%  Sample  the  spreading  eodes 
ebest  =  kron(best,wf); 
eworst  =  kron(worst,wf); 

%  Caleulate  the  signal  power 

%  The  rbw  in  the  signal  has  unit  power  multiple  effeet 
%  thus  the  power  is  just  the  power  of  the  spreading  eode 
%  All  three  spreading  eodes  have  the  same  power 
sigpower  =  ebest(l,:)*ebest(l,:)'/(samps*slen); 

%%%%%%  Normalize  the  energy  in  the  signals 
ebest  =  ebest./sqrt(sigpower); 
eworst  =  eworst./sqrt(sigpower); 

sigpowern  =  ebest(l,:)*ebest(l,:)'/(samps*slen); 

%%%%  Gale  noise  power  required  for  Eb/No 
%  noise  eoeef  =  0.5*samps*slen*(10^(eb/10))^(-l) 
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npower  =  0.5*sigpower_n*samps*slen*(10.^(eb/10)).^(-l); 
for  ii  =  1  :length(eb) 

errorsb  =  0;  %%%%  Zero  out  number  of  errors 

errorsw  =  0; 
errorsr  =  0; 

totaltested  =  0;  %%%  Zero  out  number  of  total  bits  tested 

totnoise  =  randn(l,slen*samps*nsamps);  %%%  Gen  noise  nsamps  number  of  bits  long 

nbit  =  1 ;  %%%  Counts  which  noise  bit  to  use 

srand  =  sign(randn(l,slen*samps*nsamps)); 

trand=  kron(srand,wf) ./ sqrt(sigpower) ; 

coden  = 1 ; 

%%%  generate  the  first  bit  for  each  user 

datal  =  sign(rand(numbsig,bits)-0.5); 

rbwl(:,l:slen*  samps)  =  datal(:,l)*ones(l,slen*samps); 

for  user  =  1  mumbsig 

crandl(user,:)  =  trand((coden-l)*slen*samps+l :  coden*  slen*  samps); 
coden  =  coden  + 1 ; 
end 

while  errorsr  <  endbit 

%%%%  Generate  new  random  spreading  codes  for  each  bit 
for  user  =  1  mumbsig 

crand2(user,:)  =  trand((coden-l)*slen*samps+l:coden*slen*samps); 
coden  =  coden  +  1 ; 
if  coden  >  nsamps 

srand  =  sign(randn(l,slen*samps*nsamps)); 
trand=  kron(srand,wf) ./ sqrt(sigpower) ; 
coden  = 1 ; 
end 
end 

crand  =  [crandl  crand2]; 

%%%  Generate  the  second  bit  for  each  user 

data2  =  sign(rand(numbsig,bits)-0.5); 

rbw2(:,l:  slen*  samps)  =  data2(:,l)*ones(l,slen*samps); 

data  =  [datal  data2]; 
rbw  =  [rbwl  rbw2]; 

%%  Number  of  bits  in  the  signal 
%  numbit  =  length(data(l))/slen; 

%%%  Noise 

%%  Create  AGWN  with  power  of  value  in  sqrt  (var  =  value  in  sqrt) 

%%  eb/no  ~  SNR*slen/2 
% 

%  Eb/No  =  0  dB  noise  has  15.5 
%  =  1  dB  noise  has  12.312 

%  =  2  dB  noise  has  9.775 

%  =  4  dB  noise  has  6. 169 

%  =  6  dB  noise  has  3.8941 

%  =  8  dB  noise  has  2.4564 

%  =  lOdB  noise  has  1.5499 
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ifnbit  >  nsamps 

totnoise  =  randn(l,slen*  samps  *nsamps);  %%%  Gen  noise  nsamps  number  of  bits  long 
nbit  =  1 ;  %%%  Counts  which  noise  bit  to  use 

end 

noise  =  sqrt(npower(ii))*totnoise(slen*samps*(nbit-l)+l  :slen*samps*nbit); 
nbit  =  nbit  +  1 ; 

clear  rbest  rworst  rrand; 
forj  =  Imumbsig 

rbest(j,:)  =  rbw(j,delay(j)+l:delay(j)+slen*samps).*cbest(l,:).*[cbest(j,delay(j)+l:slen*samps) 
cbest(j,l  :delay(j))];  %%%  spread  and  despread  symbols 

rworst(j,:)  =  rbw(j,delay(j)+l  :delay(j)+slen*samps).*cworst(l,:).*[cworst(j,delay(j)+l  :slen*samps) 
cworst(j,l : delay (j))];  %%%  spread  and  despread  symbols 
rrand(j,:)  = 

rbw(j,delay(j)+l:delay(j)+slen*samps).*crand(l,l:341).*crand(j,delay(j)+l:delay(j)+slen*samps);  %%% 
spread  and  despread  symbols 

end 

rbest(numbsig+l,:)  =  noise. *cbest(l,:); 
rbest(numbsig+2,:)  =  sum(rbest); 

rworst(numbsig+l,:)  =  noise. *cworst(l,:); 
rworst(numbsig+2,:)  =  sum(rworst); 

rrand(numbsig+l,:)  =  noise. *crand(l, 1:341); 
rrand(numbsig+2,:)  =  sum(rrand); 

rxbitb  =  sign(sum(rbest(numbsig+2,l:slen*samps),2)); 

rxbitw  =  sign(sum(rworst(numbsig+2,l:slen*samps),2)); 

rxbitr  =  sign(sum(rrand(numbsig+2,l:slen*samps),2)); 

totaltested  =  totaltested  +  bits; 

errorsb  =  errorsb+abs((data(l,l)-rxbitb)/2); 
errorsw  =  errorsw+abs((data(l,l)-rxbitw)/2); 
errorsr  =  errorsr+abs((data(l,l)-rxbitr)/2); 

if  errorsr  >=  endbit 
break 
end 

datal  =  data2; 
rbwl  =  rbw2; 
crandl  =  crand2; 

end 

berb(ii)  =  errorsb/totaltested; 
berw(ii)  =  errorsw/totaltested; 
berr(ii)  =  errorsr/totaltested; 

end 


A-4 


A.2  ber_shape_sync_BR51  l.m 

function  [berb,berr]  =  ber_shape_sync_BR511(eb,numbsig, samps, shape) 

% 

%  [berb  berr]  =  ber_shape_sync_BR51  l(eb,numbsig, samps, shape) 

% 

%  Experimentally  determine  the  BER  for 
%  numbsig  users  at  each  of  the  eb  values 
%  of  Eb/No  (dB).  All  signals  are  synchronous. 

% 

%  Inputs:  eb  =  Eb/No  (dB) 

%  numbsig  =  number  of  signals  (including  desired  user) 

%  samps  =  samples  per  chip  interval 
%  srand  =  random  codes 
%  shape  =  desired  shape  for  the  chip  interval 
%  1  =  rectangular 

%  2  =  half  sine 

%  3  =  raised  cosine 

%  4  =  Blackman 

%  5  =  Lanczos 

% 

%  Returns 

%  berb  =  BER  using  best  case  Gold  codes 
%  berr  =  BER  using  random  codes  defined  by  srand 
% 

% 

%  Written  by  Matthew  Glen  5  Dec  03 


%%%%  Thesis  Code 

bits  =  1 ; 
slen  =  511; 
nsamps  =  50; 
t  =  0:samps-l; 

if  shape  ==  1 
wf=  ones(l, samps); 
elseif  shape  ==  2 

wf  =  sin(pi*t/(length(t)-l)); 

elseif  shape  ==  3 

wf  =  sqrt(2/3)*(l  -  cos(2*pi*t/(length(t)-l))); 
elseif  shape  ==  4 

wf  =  0.42  -  0.5*cos(2*pi*t/(length(t)-l)); 

elseif  shape  ==  5 
t2  =  t+10; 

wf  =  (sin(2*pi*t2)./(2*pi*t2)).^2; 
end 

%%%%%  Create  m-seq  needed  for  Gold  codes 

gpolyl  =  [1  0000  1  000  1];  %%%  1021  octal  polynomial 

gpoly2  =  [1  1001  10001];  %%%  1461  octal  polynomial 
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iiiit=  [111111111];  %%%  initial  register  state 

mseq  =  [zeros(2,5 11)];  %%%  create  mseq  vector  of  zeros 

mseq(l,:)  =  mSeqGen(init,gpolyl);  %%%  generate  the  first  mseq 
mseq(2,:)  =  mSeqGen(init,gpoly2);  %%%  generate  the  2nd  mseq 


gold  =  zeros(slen+2,slen); 

gold  =  gold_gen(mseq(l,:),  mseq(2,:)).*(-2)  +1;  %%  Generate  the  Gold  Code  Family 
[row  col]  =  size(gold); 


%%  Run  test  until  endbit  number  of  errors  is  reached 
endbit  =  300; 

%  Define  the  best  and  worst  combinations  of  51 1-length  gold  codes 
bestorder  =  [1  3:51 1];  %  best  codes 

%  worstorder  =  [2  8  11  12  20  23  25  27  29  31  32  34  36  37  41  42  43  45  46  47  48  49  51 
%  52  57  60  61  62  68  72  79  80  81  84  94  96  98  100  101];  %  worst  codes 

%%%  Spreading  codes 

best(l:numbsig,:)  =  gold(bestorder(l:numbsig),:); 


%  Sample  the  spreading  codes 
cbest  =  kron(best,wf); 

%  Make  2  bits  worth  of  the  spreading  code 
%  cbest  =  [cbest  1  cbest  1]; 

%  cworst  =  [cworstl  cworstl]; 

%  crand  =  [crandl  crandl]; 


%  Calculate  the  signal  power 

%  The  rbw  in  the  signal  has  unit  power  multiple  effect 
%  thus  the  power  is  just  the  power  of  the  spreading  code 
%  All  three  spreading  codes  have  the  same  power 
sigpower  =  cbest(l,:)*cbest(l,:)'/(samps*slen); 

%%%%%%  Normalize  the  energy  in  the  signals 
cbest  =  cbest./sqrt(sigpower); 

%%  Calculate  new  Normalized  signal  power 

%%  This  value  remains  the  same  as  sigpower  for  rectangular  shapes 

sigpowern  =  cbest(l,:)*cbest(l,:)'/(samps*slen); 

%%%%  Calc  noise  power  required  for  Eb/No 
%  noise  coeef  =  0.5*samps*slen*(10^(eb/10))''(-l) 
npower  =  0.5*sigpower_n*samps*slen*(10.''(eb/10))A(-l); 

%%  Creates  a  vector  of  zeros  to  be  fdled  with  the  random  spreading  codes 
crand  =  zeros(9,slen*samps); 


for  ii  =  1  :length(eb) 

errorsb  =  0;  %%%%  Zero  out  number  of  errors 
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errorsr  =  0; 

totaltested  =  0;  %%%  Zero  out  number  of  total  bits  tested 

%%%  Generate  nsamps  number  of  bits  worth  of  noise 
totnoise  =  randn(l,slen*samps*nsamps); 

nbit  =  1 ;  %%%  Counts  whieh  noise  bit  to  use 

%%%  Generate  a  string  of  random  I's  and  -I's  to  be  used  as  the  random 

%%%  spreading  codes 

srand  =  sign(randn(l,slen*samps*nsamps)); 

%%  Track  the  number  of  bits  worth  of  random  spreading  code  that  has 
%%  been  used 
coden  =  1 ; 

while  errorsb  <  endbit  %%  Run  until  'endbit'  number  of  errors 

%%%%  Generate  new  random  spreading  codes  for  each  bit 

%%  Create  the  random  spreading  codes  for  each  transmitting  user 
for  user  =  1  mumbsig 

crand(user,:)  =  kron(srand((coden-l)*slen+l  :coden*slen),wf)./sqrt(sigpower); 
coden  =  coden  + 1 ; 

%%  Generate  new  random  1  's  and  -  Ts  when  the  previous  srand  is 
%%  all  used  up 
if  coden  >  nsamps 

srand  =  sign(randn(l,slen*samps*nsamps)); 
coden  = 1 ; 
end 
end 

%%  Generate  the  current  data  bit  values  for  each  user 
data  =  sign(rand(numbsig,bits)-0.5); 

%%  Sample  the  current  data  bit  values  for  each  user 
rbw(:,l:slen*  samps)  =  data(:,l)*ones(l,slen*samps); 

%%%  Noise 

%%  Create  AGWN  with  power  of  value  in  sqrt  (var  =  value  in  sqrt) 

%%  eb/no  ~  SNR*slen/2 
% 

%  Eb/No  =  0  dB  noise  has  15.5  variance 
%  =  1  dB  noise  has  12.312  variance 

%  =  2  dB  noise  has  9.775  variance 

%  =  4  dB  noise  has  6. 169  variance 

%  =  6  dB  noise  has  3.8941  variance 

%  =  8  dB  noise  has  2.4564  variance 

%  =  lOdB  noise  has  1.5499  variance 

%%  Group  a  symbol's  worth  of  noise  samples  to  be  added  to  received 
%%  signal  from  nsamps  long  noise  string  'totnoise'.  Generate  new 
%%  long  noise  string  if  nsamps  number  of  bits  have  already  been 
%%  tested, 
if  nbit  >  nsamps 

totnoise  =  randn(l,slen*  samps  *nsamps);  %%%  Gen  noise  nsamps  number  of  bits  long 
nbit  =  1 ;  %%%  Counts  which  noise  bit  to  use 

end 
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noise  =  sqrt(npower(ii))*totnoise(slen*samps*(nbit-l)+l  :slen*samps*nbit); 
nbit  =  nbit  +  1 ; 

clear  rbest  rrand; 
forj  =  l:numbsig 

rbest(j,:)  =  rbw(j,:).*cbest(l,:).*cbest(j,:);  %%%  spread  and  despread  symbols 

rrand(j,:)  =  rbw(j,:).*crand(l,:).*crand(j,:);  %%%  spread  and  despread  symbols 

end 

%%  Current  received  signal  for  Best-Case  Gold  Coded  signals 
rbest(numbsig+l,:)  =  noise. *cbest(l,:); 
rbest(numbsig+2,:)  =  sum(rbest); 

%%  Current  received  signal  for  randomly  Coded  signals 
rrand(numbsig+l,:)  =  noise. *crand(l,:); 
rrand(numbsig+2,:)  =  sum(rrand); 

%%  Estimate  the  current  "best-case"  received  bit 
rxbitb  =  makeint(sum(rbest(numbsig+2,l:slen*samps),2)); 

%%  Estimate  the  current  "random"  received  bit 

rxbitr  =  makeint(sum(rrand(numbsig+2,l:slen*samps),2)); 

totaltested  =  totaltested  +  bits; 

%%  Determine  if  estimated  bits  are  in  error  and  increment  total 

%%  number  of  errors 

errorsb  =  errorsb+abs((data(l,l)-rxbitb)/2); 

errorsr  =  errorsr-l-abs((data(l,l)-rxbitr)/2); 

end 

berb(ii)  =  errorsb/totaltested; 
berr(ii)  =  errorsr/totaltested; 
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